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CHAPTER  I 


INTRODUCTION 

The  basic  subject  of  this  thesis  is  setting  optimum 
safety  stock  levels  in  multiple- item  vendor  inventories  under 
system  constraints.  An  inventory  is  essentially  an  idle  re¬ 
source  which  is  being  temporarily  stored  for  use  at  some 
future  time  (16:1).  The  primary  purpose  for  storing  resources 
to  decouple  the  otherwise  dependent  functions  in  the 
production-distribution-consumption  chain  (2:389)  by  separa¬ 
ting  the  supply  and  demand  processes  (16:1).  Vendor  (pur¬ 
chasing)  inventories  are  those  which  are  procured  from  outside 
the  organization  as  opposed  to  production  inventories  which 
are  produced  and  stored  by  the  eventual  user  (2:390).  The 
function,  then,  of  vendor  inventories  is  to  separate  the 
supply  process  involved  in  procuring  the  resource  from  the 
demand  process  of  consuming  the  resource  or  redistributing 
(selling)  the  resource  outside  the  organization.  Vendor  in¬ 
ventories  are  typically  multiple-item  inventories  in  that 
most  organizations  store  several  different  resources  rather 
than  a  single  item  (2:391). 

Because  maintaining  any  type  of  inventory  incurs 
expenses,  the  benefits  of  holding  the  inventory  must  exceed 
the  holding  costs  (16:6).  The  process  of  balancing  the  costs 
incurred  with  the  costs  avoided  by  holding  stocks  is  inventory 
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control  (16:7).  As  Hadley  and  Whitin  (9:1)  state,  the  two 
fundamental  issues  in  controlling  any  inventory  are  when  to 
order  and  how  much  to  order.  In  addition,  when  the  demand 
for  an  item  is  uncertain,  the  level  of  safety  or  buffer  stock 
carried  to  meet  unpredicted  demand  must  be  considered  when 
making  these  two  decisions  (9:161).  The  specific  subject  of 
this  thesis  is  the  problem  of  determining  the  optimum  item 
safety  stock  levels  in  the  Defense  Electronics  Supply  Center 
(DESC)  inventory  subject  to  some  constraints.  The  type  of 
constraints  considered  are  those  that  place  limits  on  the 
combined  safety  stock  level  of  all  the  items  in  the  DESC  in¬ 
ventory  . 

The  background  for  the  research  problem  is  presented 
in  the  next  section,  which  defines  the  operating  characteris¬ 
tics  of  the  DESC  inventory  and  describes  the  organizational 
setting  for  the  problem. 


Background 

The  DESC  Inventory 

DESC  manages  and  controls  the  inventory  of  non¬ 
reparable  electronics  parts  and  components  for  the  U.S.  De¬ 
partment  of  Defense  (DoD) .  DESC  procures  this  inventory  from 
numerous  sources  and  then  supplies  the  various  DoD  components 
with  items  as  requested  (23).  Thus,  the  DESC  inventory  is  a 
vendor  inventory  in  which  items  are  replenished  instantaneously 
when  an  order  arrives  from  a  supplier.  Vendor  inventories  are 
normally  differentiated  from  production  inventories  which  are 
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replenished  gradually  as  the  items  are  produced  (2:390). 
Another  common  distinction  made  between  types  of  inventories 
is  whether  they  are  static  or  dynamic.  Static  inventories 
are  those  for  which  a  single  procurement  is  made  to  satisfy 
the  demand  in  a  single  finite  time  period,  while  a  dynamic 
inventory  has  an  infinite  planning  horizon  and  the  inventory 
is  repeatedly  restocked  (24:6).  DESC  restocks  items  as  re¬ 
quired;  therefore,  the  DESC  inventory  can  be  classified  as 
a  dynamic  vendor  inventory. 

Inventories  may  also  differ  in  how  often  the  stock 
level  is  checked.  Stocks  can  be  reviewed  either  continuously 
with  fixed-quantity  orders  placed  at  any  time  the  reorder 
point  is  reached,  or  they  can  be  reviewed  periodically  to 
determine  when  and  how  much  to  order  (16:9).  Although  DESC 
reviews  stock  levels  at  three-day  intervals,  the  inventory  is 
controlled  as  if  it  were  a  continuous  review  inventory  (1). 
This  is  reasonable  because  continuous  review  can  be  assumed 
for  practical  purposes  if  the  time  interval  between  reviews 
is  negligible  compared  to  the  interval  between  successive 
stock  depletions  (16:40).  Figure  1  shows  the  basic  behavior 
for  a  single  item  in  a  continuous  review,  dynamic,  vendor 
inventory.  The  inventory  level  decreases  over  time  at  the 
rate  of  demand,  D,  until  an  order  is  received.  Then  the  stock 
is  instantaneously  replenished  by  the  standard  quantity,  Q, 
which  was  ordered  when  the  inventory  level  reached  the  re¬ 
order  point,  R. 

The  DESC  inventory  can  be  further  classified  as  one 
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which  allows  backorders  as  opposed  to  the  lost  sales  type. 
Inventories  which  allow  backorders  eventually  satisfy  all 
shortages  as  stocks  are  replenished.  If  backorders  are  not 
allowed,  any  demand  that  cannot  be  satisfied  from  current 
stocks  represents  lost  sales  (2:391).  Figure  2  depicts  the 
behavior  of  an  inventory  item  resulting  from  allowing  back¬ 
orders  in  the  model  shown  in  Figure  1.  In  the  model  with 
backorders,  the  order  quantity  is  used  both  to  fill  backorders 
and  to  replenish  the  inventory.  The  advantage  of  allowing 
backorders  is  that  fewer  orders  are  placed  and  inventory 
levels  are  lower.  Consequently,  both  the  costs  of  ordering 
stocks  and  holding  stocks  may  be  reduced  (2:401). 

The  final  characteristic  to  be  considered  in  defining 
the  DESC  inventory  is  whether  demands  and  leadtimes  are  deter¬ 
ministic  or  stochastic.  The  models  shown  in  both  Figures  1 
and  2  assume  that  the  demand  for  an  item  and  the  delay,  or 
leadtime,  between  placing  an  order  and  receiving  the  order 
are  known  with  certainty,  and  that  demand  is  uniform  and  con¬ 
tinuous  with  respect  to  time.  However,  the  demand  or  leadtime 
may  be  variable  and  unknown  and  would,  therefore,  be  stochas¬ 
tic  rather  than  deterministic.  When  the  demand  or  leadtime 
are  stochastic,  safety  stocks  can  be  maintained  to  absorb 
some  amount  of  variation  in  the  demand  or  leadtime  (2:390). 

In  the  DESC  system,  safety  stocks  must  be  carried  because  both 
the  demand  and  leadtime  are  stochastic.  Figure  3  is  a  model 
of  this  situation.  Although  backorders  are  allowed,  safety 
stocks  are  still  necessary  to  insure  excessive  backorders  do 


Time 


II 


Figure  2 


Behavior  of  a  Continuous  Review, 
Dynamic,  Vendor  Inventory  with 
Backorders 


Figure  3 


Inventory  with  Stochastic 
Demand  and  Leadtimes 


not  accumulate. 

Although  the  DESC  inventory  system  could  be  defined 
further,  the  characteristics  described  are  sufficient  to 
understand  the  environment  in  which  the  optimum  safety  stock 
levels  must  be  determined.  Backorders  are  allowed  for  each 
of  the  numerous  different  items  in  the  inventory.  For  each 
item,  safety  stock  can  be  carried  to  either  avoid  shortages 
or  to  limit  the  number  of  backorders.  Thus,  if  high  levels 
of  safety  stock  are  carried  for  each  item,  the  total  number 
of  backorders  is  reduced  and  the  overall  system  performance 
is  improved  in  terms  of  being  able  to  supply  required  items 
to  the  DoD  components.  At  the  same  time,  higher  safety  stock 
levels  require  a  larger  inventory  investment.  The  problem, 
then,  is  to  determine  the  optimum  safety  stock  level  for  each 
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inventory  item  based  on  the  level  of  system  performance  which 
can  be  achieved  with  the  available  funds.  The  organizational 
setting  of  this  problem  is  briefly  described  in  the  next  sec¬ 
tion. 

DESC  Organization 

DESC  is  a  component  of  the  Defense  Logistics  Agency 
(DLA) ,  a  DoD  organization  responsible  for  providing  services 
and  supplies  which  are  used  in  common  by  all  the  Military 
Services.  DLA  was  established  in  1962  to  consolidate  several 
DoD  supply  agencies  which  had  been  managed  by  separate  Mili¬ 
tary  Services.  Today,  DESC  is  one  of  six  Defense  Supply 
Centers  under  the  direction  of  DLA  (25:1).  DESC  is  the  DoD 
inventory  control  point  for  non- reparable  electronics  parts 
and  components.  DESC  is  responsible  for  the  supply  management 
of  electronics  components  used  by  the  Military  Services,  in¬ 
cluding  consolidating  requirements,  procurement,  storage, 
distribution,  and  financial  accounting.  In  addition,  DESC 
provides  supply  support  to  all  Federal  civil  agencies  (27:11-2). 

DESC  is  presently  the  inventory  control  point  for 
over  750,000  items  which  represent  an  annual  inventory  invest¬ 
ment  of  over  $450  million.  This  inventory  supplies  some 
22,000  military  and  civil  agency  customers  (23).  In  control¬ 
ling  this  inventory,  DESC  is  responsible  for  implementing  DoD 
policy  to  determine  order  quantities,  reorder  points,  and 
other  inventory  operating  rules  for  each  item.  Because  the 
supply  and  demand  processes  differ  greatly  amongst  the  various 
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items,  DESC  uses  different  control  procedures  for  different 
categories  of  items  based  on  such  factors  as  the  demand  rate, 
leadtime,  and  item  cost.  Of  specific  interest  here  is  that 
some  items  are  controlled  using  fixed  safety  stock  levels, 
while  others  have  variable  safety  levels  which  are  changed 
periodically  based  on  updated  information  (1).  This  study  is 
concerned  only  with  the  total  backorders  and  the  associated 
funding  requirement  for  those  items  which  DESC  controls  by 
using  a  variable  safety  level  operating  rule. 

DLA/DESC  Variable  Safety  Level 

The  DoD  policy  and  the  methods  used  by  DLA  and  DESC 
to  determine  safety  stock  levels  have  evolved  and  changed  over 
the  years.  This  section  presents  the  background,  development, 
and  application  of  the  model  presently  used  by  DESC  and  DLA 
to  determine  variable  safety  levels. 

DoD  Policy 

Safety  level  is  defined  by  DoD  as  "the  quantity  of 
material  which  is  required  to  be  on  hand  to  permit  continued 
operation  in  the  event  of  minor  interruption  of  normal  replen¬ 
ishment  or  unpredictable  fluctuation  in  demand  [28:2]." 
According  to  a  DLA  inventory  management  specialist  (25:1), 
the  only  DoD  guidance  on  the  subject  of  safety  levels  prior 
to  1970  was  DoD  Instruction  4140.11,  which  simply  stated  that 
variable  safety  levels  were  preferable  to  fixed  levels.  Each 
DoD  component  was  left  with  the  responsibility  to  establish 


an  appropriate  policy  and  develop  mathematical  models  to  deter¬ 
mine  inventory  safety  levels.  As  a  result,  DoD  supply  centers 
were  operating  under  a.  wide  range  of  safety  level  policies  of 
varying  effectiveness  and  cost. 

In  1970  the  DoD  published  its  policy  for  determining 
order  quantities  and  safety  stock  levels  for  non- reparable 
items  in  DoD  Instruction  4140.39.  This  policy  specifies  both 
a  system  objective  function  which  should  be  optimized  and  a 
system  constraint  which  must  be  met.  The  objective  function 
is  the  total  variable  cost  of  ordering  and  holding  inventory 
at  the  inventory  control  point.  This  cost  should  be  minimized 
subject  to  a  system  constraint  on  the  number  of  requisitions 
short  (backorders)  considering  the  average  number  of  days 
delay  and  the  essentiality  of  each  item.  Thus,  the  objective 
of  this  policy  is: 

To  minimize  the  total  of  variable  order  and  holding 
costs  subject  to  a  constraint  on  time-weighted, 
essentiality-weighted  requisitions  short  [28:2]. 

DoD  Instruction  4140.39  also  prescribed  the  general 
model  to  be  used  in  implementing  the  policy  statement.  The 
annual  variable  order  cost  for  each  item  was  given  as: 


where 

A  *  order  cost 

‘f  V\ 

D^  *  mean  annual  demand  for  the  il  item 
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where 

I  =  holding  cost  rate 
=  cost  of  item 

The  expected  essentiality-weighted  requisitions  short  (in  terms 
of  orders  rather  than  units)  at  a  random  point  in  time  was 
expressed  as: 
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RS  =  Z  ,*-4-/“  (x-R.)[F(x+Q.  ; L)  -  F(x;L)]dx 
i= 1  Ki  1  1 

where 

»  item  essentiality 

=  average  units  demanded  per  requisition 

F(x+Q^;L)  =  probability  that  the  number  of  units  demanded 
in  leadtime  L  is  <  x  +  Q- 

F(x;L)  *  probability  that  the  number  of  units  demanded 
in  leadtime  L  is  <_  x. 

Finally,  the  inventory  total  variable  cost  to  be  minimized  was 
defined  as : 

TVC  =  OC  +  HC  +  ARS 

where  XRS  is  the  implied  cost  of  time-weighted  shortages. 
Minimizing  this  total  variable  cost  is  equivalent  to  using  the 
method  of  Lagrange  multipliers  to  minimize  the  sum  of  the  order¬ 
ing  and  holding  cost  subject  to  a  constraint  on  the  expected 
es sentiality- weighted  requisitions  short  (11:552). 

The  expected  backorder  model  used  to  determine  the 
implied  shortage  cost  was  derived  from  a  model  developed  by 
Hadley  and  Whitin  in  1965  (25:3).  The  implementation  of  the 
model,  including  determination  of  the  shortage  parameter  (X), 
the  choice  of  leadtime  demand  distributions,  and  the  develop¬ 
ment  of  explicit  formulas  was  left  to  the  individual  DoD  com¬ 
ponents.  The  instruction  noted  that  if  the  shortage  parameter 
is  specified  in  terms  of  a  performance  level,  this  will  affect 
the  required  funding  level  and  vice  versa.  Further,  it  pointed 
out  that  the  shortage  parameter  used  by  the  DoD  to  compute  the 
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budget  would  be  a  function  of  performance  goals  while  the  para¬ 
meter  used  in  day-to-day  operations  might  be  a  function  of 
funding  levels  or  other  management  decisions  (28£ncl  2,  p.  3). 


The  DLA  Safety  Level  Model 

In  implementing  the  DoD  policy,  DLA  adopted  a  modified 
Wilson  EOQ  model  to  determine  order  quantities  and  a  model 
developed  by  Presutti  and  Trepp  to  determine  safety  stock 
levels  (6:A-1).  The  model  chosen  to  determine  item  safety 
levels  uses  the  number  of  backorders  (a  performance  measure) 
as  the  system  constraint,  and  assumes  that  demands  are  norm¬ 
ally  distributed  using  the  following  as  an  approximation  to 
the  normal  probability  density  function  (18  :244) : 

F(x)  =  exp(-/Ii^|);  for  all  x 

Using  a  technique  developed  by  Hadley  and  Whitin  (9:178), 
Presutti  and  Trepp  show  (18  :246)  that,  assuming  the  above  func¬ 
tion  describes  the  demand  distribution,  the  expected  number 
of  units  in  a  backorder  status  at  a  random  point  in  time  is: 

2 

Bt  *  (l-exp(-/7  Q/a))  exp  (-/Ik) 

Presutti  and  Trepp  then  use  this  expression  to  develop  four 
inventory  models  which  differ  in  their  treatment  of  backorder 
penalties  and  holding  costs.  DLA  chose  to  adopt  the  model  in 
which  the  backorder  penalty  is  time -weighted  and  the  holding 
cost  is  applied  to  the  inventory  position  (on-hand  plus  due- ins) 
rather  than  the  on-hand  inventory  (6:A-2)  used  in  the  DoD  model. 
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The  basic  DLA  model,  using  DLA  notation  is  (6:A-2): 

N  P . AD .  N  Q . 

minimize  Z  — —  +  E  a.C.(u.  +  k.o-  +  J1) 

i=l  yi  i=l  1  1  1  11  1 

subject  to 

N  .  Z . o • 2  Q. 

z  -Hr  Cl-exp(-/J  -±))  exp {-/l  k  )  <  0 

i-1  z  i^i  ai  1  " 

where 

k^  =  safety  level  factor 

=  Economic  Order  Quantity 

=  leadtime  demand  or  expected  due- ins  (3:187) 
a^  =  holding  cost  rate 
=  unit  price 

=  standard  deviation  of  leadtime  demand 
=  essentiality  factor 
=  average  requisition  size 
AD^  =  annual  demand  in  units 
=  procurement  (order)  cost 
8  =  number  of  backorders  at  a  random  point  in  time 

Consistent  with  the  DoD  model,  this  model  minimizes  total  order 
ing  and  holding  costs  subject  to  a  constraint  on  the  number  of 
backorders.  The  Method  of  Lagrange  then  yields  the  following 
formula  to  compute  item  safety  level  factors: 


where 


k. 


i 


In  [- 


St/2  QjajCj 


.  5(-X)Ziai(l-exp(-/2 
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-  X 


N  aiaiCi 
i  =  l  ^6 


DLA  uses  an  approximation  for  the  standard  deviation  of  lead- 

time  demand,  which  is  1.25  times  the  mean  absolute  deviation 

of  demand  forecast  errors  over  the  leadtime  (o^  =  1.25  MADLT) 

and  a  constant  holding  rate  (a^  =  a) .  Using  these  values  and 

substituting  the  expression  for  X  into  k^  yields: 

t  2.56  SiQiCi6 

ki  =  '  7f  ln[Z~MAl)TTiTETi  m'DLTinTr‘iTpTr72'0i71. 2'$'  MaCLT.)!1 

The  variable  safety  level  for  an  item  is  then  calculated  as 


VSL .  =  1.25  k.  MADLT. 
l  ii 

which  is  the  safety  level  factor  times  the  approximation  for 
the  standard  deviation  of  the  leadtime  demand. 


Model  Application 

DESC  has  used  the  above  model  to  calculate  variable 
safety  levels  for  about  150,000  inventory  items  each  quarter 
since  January  1976.  Other  items,  such  as  new  items  with  no 
demand  history  or  items  requisitioned  one  time,  are  managed 
using  various  other  operating  rules.  From  January  1976  through 
December  1980,  the  average  investment  in  variable  safety  level 
stock  was  about  $53  million,  and  the  standard  error  from  the 
desired  funding  level  was  about  $5.9  million.  During  this 
period,  the  method  used  to  include  the  desired  funding  level 
as  a  constraint  on  safety  stock  levels  was  to  use  experience 
to  select  a  backorder  constraint  which  would  result  in  the 
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desired  safety  level  dollars  after  the  item  safety  levels 
were  calculated  according  to  the  DLA  model.  Thus,  the  objec¬ 
tive  in  selecting  a  backorder  constraint  as  an  input  to  the 
DLA  model  has  been  to  choose  a  value  which  sets  item  safety 
levels  such  that  the  total  required  safety  stock  investment 
equals  the  desired  funding  level  (1)  . 

The  selected  backorder  constraint  was  then  input 
to  a  computer  program  which  calculates  the  item  variable 
safety  levels  based  on  the  DLA  model.  This  same  program  also 
accomplishes  the  inventory  quarterly  update,  calculates 
numerous  other  item  values,  and  after  running  about  18  hours, 
produces  a  new  inventory  policy  to  be  used  in  the  upcoming 
quarter  for  each  item.  In  addition  to  computing  the  safety 
level  dollars  required  for  the  input  backorder  constraint, 
the  program  also  computes  five  other  funding  levels  based  on 
different  (selected]  backorder  constraints.  The  performance 
values  (backorders)  are  then  plotted  against  the  funding 
levels  (safety  level  dollars)  to  derive  a  curve  which  shows 
the  correct  backorder  constraint  which  should  have  been  used 
to  achieve  the  funding  target.  However,  to  use  the  correct 
constraint  for  the  upcoming  quarter  would  require  rerunning 
the  complete  quarterly  update  program  which  would,  of  course, 
be  extremely  costly  due  to  the  amount  of  computer  time  re¬ 
quired.  Therefore,  the  post- analysis  to  determine  the  correct 
backorder  constraint  has  been  of  little  practical  use  (1)  . 
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Problem  Statement 


At  the  present  time,  the  model  used  by  DESC ,  and  all 
other  DLA  supply  centers,  to  determine  variable  safety  levels 
does  not  explicitly  consider  any  constraint  on  funds  available 
for  the  inventory  safety  stoclc.  The  funding  constraint  is 
implicitly  included  in  the  model  by  selecting  a  backorder 
level  which  will  maximize  system  performance  while  satisfying 
the  funding  constraint.  Using  experience  to  select  the  back¬ 
order  constraint  has  produced  relatively  large  errors  between 
the  target  funding  level  and  the  required  safety  level  dollars. 
This  causes  the  individual  item  variable  safety  levels  for  a 
particular  quarter  to  be  set  at  a  level  which  is  significantly 
different  from  the  optimum  level  considering  the  funds  avail¬ 
able.  Thus,  the  basic  problem  of  the  current  model  is  that 
it  is  not  effective  in  establishing  safety  levels  that  are 
consistent  with  funding  targets.  DESC  is  interested  in  re¬ 
solving  this  problem  to  gain  better  control  of  the  system  and 
improve  service  (1) . 

Research  Objectives 

Although  DESC  inventory  managers  are  concerned  about 
maintaining  certain  performance  levels,  the  primary  constraint 
in  establishing  item  safety  levels  is  the  funds  available  for 
safety  stock.  Therefore,  the  overall  objective  of  the  research 
was  to  develop  an  efficient  and  accurate  method  to  determine 
item  safety  levels  subject  to  a  budget  constraint.  However, 


since  an  accurate  prediction  of  system  performance  is  impor¬ 
tant  and  can  be  used  to  demonstrate  the  effects  of  different 
funding  levels,  another  goal  was  to  determine  the  relationship 
between  performance  levels  and  funding  requirements.  In  fact, 
if  this  relationship  is  established  prior  to  running  the 
quarterly  update,  then  it  will  make  no  difference  whether  the 
system  constraint  used  in  the  model  is  in  terms  of  performance 
level  or  funding  level.  The  primary  objective  of  the  research, 
then,  was  to  develop  a  method  to  define  the  relationship,  at  a 
specific  point  in  time,  between  the  expected  number  of  total 
backorders  and  the  required  funding  level.  This  information 
can  then  be  used  to  set  item  safety  levels  based  on  the  level 
of  system  performance  which  can  be  achieved  with  the  available 
funds . 


Summary 

This  chapter  introduced  the  research  problem.  After 
some  preliminary  comments  about  inventories  in  general,  the 
DESC  inventory  system  was  defined  as  a  continuous  review, 
dynamic,  vendor  inventory  system  which  allows  backorders  and 
is  subject  to  stochastic  demands  and  leadtimes.  As  a  part  of 
the  Defense  Logistics  Agency  organization,  DESC  uses  the  DLA 
prescribed  model  to  determine  item  variable  safety  levels. 

The  objective  of  this  model  is  to  minimize  the  total  system 
costs  of  ordering  and  holding  the  inventory  subject  to  a  con¬ 
straint  on  the  total  number  of  backorders  allowed.  The  prob¬ 
lem  being  experienced  with  this  model  is  that  it  does  not 


consider  the  amount  of  funds  available  for  safety  stock, 
which  causes  large  errors  between  the  required  safety  level 
investment  and  the  funding  target. 
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CHAPTER  2 


PROBLEM  DISCUSSION 

The  research  problem,  as  introduced  in  Chapter  1,  is 
to  determine  the  relationship  between  total  inventory  back¬ 
orders  and  total  safety  stock  investment  when  the  variable 
safety  stock  levels  are  set  using  the  current  DLA/DESC  model. 
This  chapter  examines  the  problem,  primarily  from  a  theoreti¬ 
cal  viewpoint,  to  establish  the  basis  for  a  solution  method. 

The  purpose  is  to  clearly  define  the  research  problem  within 
the  framework  of  generally  accepted  inventory  principles,  the 
objectives  of  the  organization,  and  the  alternatives  available 
to  solve  or  alleviate  the  problem.  The  first  section  presents 
a  background  for  the  problem  by  identifying  the  inventory 
characteristics  which  lead  to  multiple  management  objectives, 
the  nature  of  the  relationships  that  exist  between  the  objec¬ 
tives,  and  the  performance  measures  used  to  evaluate  multiple- 
item  inventory  objectives.  Next,  the  theory  supporting  the 
DLA  model  is  examined  to  identify  the  strengths  and  weaknesses 
of  the  model  and  the  source  of  the  research  problem.  The  third 
section  discusses  the  problem  of  finding  the  costs  for  main¬ 
taining  different  levels  of  backorders  and  presents  the  basic 
alternatives  available  to  overcome  the  problem.  Finally,  the 
last  section  briefly  presents  some  alternative  models  for 
calculating  safety  stock  levels  which  might  provide  a  solution 
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to  the  problem  if  DESC  were  willing  to  adopt  a  new  model. 

Safety  Stock  Management  Objectives 

As  suggested  in  Chapter  1,  inventory  control  is  con¬ 
cerned  with  balancing  the  costs  and  the  benefits  of  holding 
stocks.  This  balance  is  directly  affected  by  the  level  of 
stock  held  by  the  organization  (14:111).  The  effects  of  low 
stock  levels  are  that  customer  demands  often  cannot  be  satis¬ 
fied,  the  number  of  stock  replenishment  orders  is  high,  many 
orders  must  be  expedited,  and  there  are  few  opportunities  for 
quantity  discounts.  Thus,  low  stock  levels  cause  low  service 
levels  and  high  ordering  costs.  High  stock  levels  will  in¬ 
crease  the  service  level  and  decrease  ordering  costs.  However, 
high  stock  levels  also  have  some  undesirable  effects  such  as 
high  storage  costs,  high  capital  investment  and  increased 
risk  of  obsolescence.  Thus,  high  stock  levels  increase  hold¬ 
ing  costs.  The  relationships  between  these  variables  are  that 
when  stock  levels  and  service  levels  increase,  holding  costs 
increase  and  ordering  costs  decrease. 

There  is  also  some  inventory  level  at  which  the 
total  costs  of  ordering  and  holding  the  inventory  are  mini¬ 
mized  (9:35).  If  the  only  objective  of  inventory  management 
were  to  minimize  the  ordering  and  holding  costs  and  these 
costs  could  be  quantified,  the  appropriate  inventory  level 
could  be  found  by  minimizing  the  sum  of  the  ordering  cost 
function  and  the  holding  cost  function.  This  inventory  level 
would,  in  turn,  determine  the  level  of  service.  The  problem 
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with  this  approach  is  that  not  all  ordering  and  holding  costs 
can  be  quantified  and  the  costs  of  inventory  shortages  must 
be  considered.  Therefore,  most  organizations  have  at  least 
two  other  inventory  management  objectives  in  addition  to  mini¬ 
mizing  the  quantifiable  ordering  and  holding  costs  (22:39). 

One  of  the  other  objectives  is  to  improve  the  level  of  service 
or  provide  a  minimum  level  of  service.  The  second  objective 
is  to  reduce  the  size  of  the  inventory  or  to  limit  the  size 
to  some  desired  level.  One  reason  for  these  two  additional 
objectives  is  that  the  associated  costs  are  difficult,  if  not 
impossible,  to  measure  directly  (16:22).  Service  level  is 
actually  a  surrogate  measure  for  shortage  costs,  while  inven¬ 
tory  size  can  be  considered  as  a  substitute  for  the  opportu¬ 
nity  cost  of  the  inventory  investment,  which  is  one  of  the 
holding  costs.  These  costs  are  always  difficult  to  measure, 
and  in  the  case  of  non-profit- seeking  organizations  such  as 
the  DoD,  it  is  often  not  possible  to  quantify  shortage  costs 
and  opportunity  costs  (31:  Chap.  7).  Thus,  the  basic  objec¬ 
tives  of  inventory  management  are  to  minimize  the  quantifiable 
ordering  and  holding  costs,  provide  some  level  of  service, 
and  limit  the  size  of  the  inventory. 

One  of  the  difficulties  in  inventory  management  is 
caused  by  the  fact  that  most  inventories  are  multiple- item 
and  there  are  interactions  among  the  items.  There  are  many 
types  of  possible  interactions  between  the  items.  For  example, 
the  items  may  be  substitutes  for  each  other;  ordering  costs 
might  be  reduced  by  ordering  several  different  items  from  the 
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same  supplier;  or  the  items  may  be  competing  for  restricted 
resources  such  as  warehouse  space  or  inventory  investment 
dollars  (9:S4).  Interaction  between  the  many  inventory  items 
requires  the  use  of  aggregate  performance  indices  and  more 
complex  models.  Aggregate  indices  of  costs,  service  levels, 
and  inventory  size  are  necessary  to  measure  the  effectiveness 
of  meeting  the  objectives  for  the  total  inventory  (29;8). 

More  complex  models  are  required  to  reflect  the  interactions 
among  the  many  different  items  (16:112). 

Another  management  problem  is  selecting  the  appro¬ 
priate  indices  and  establishing  specific  inventory  objectives 
in  terms  of  the  selected  indices.  There  are  numerous  indices 
in  common  use  that  provide  aggregate  measures  of  service  level 
or  inventory  size.  Some  of  these  are  (29:8): 

1.  Value  of  inventory  on  hand 

2.  Value  of  inventory  on  order 

3.  Value  of  inventory  on  hand  and  on  order 

4.  Number  of  orders  placed 

5.  Number  of  backorders 

After  selecting  indices  which  are  appropriate  for  the  parti¬ 
cular  organization,  the  management  objectives  must  be  speci¬ 
fied  in  terms  of  the  indices.  However,  the  ordering  and  hold¬ 
ing  costs,  service  level,  and  inventory  size  are  not  indepen¬ 
dent  of  each  other.  Because  these  objectives  are  interrelated, 
they  must  be  considered  simultaneously  and  tradeoffs  made  to 
achieve  a  balance  which  is  considered  optimal  for  the  parti¬ 
cular  organization  in  question. 
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The  selected  balance  between  the  desired  objectives 
may  not  be  achieved  because  of  the  stochastic  processes  inher¬ 
ent  in  many  inventory  systems.  The  distinguishing  characteris¬ 
tic  of  stochastic  inventories  is  that  stocks  are  not  ordered 
to  meet  a  particular  known  demand  (16:56).  Instead,  stocks 
are  ordered  based  on  the  current  inventory  status  and  the 
knowledge  of  past  demands  and  leadtimes  with  uncertain  future 
demands  and  leadtimes.  Because  the  demand  and  leadtime  for 
single  items  is  uncertain,  the  inventory  size,  service  level, 
and  costs  are  probabilistic  and  can  be  viewed  as  random  vari¬ 
ables.  Indices  which  are  based  on  many  items  with  probabilis¬ 
tic  demand  or  leadtime  must  be  described  stochastically  with 
the  leadtime  and  demand  pattern  of  the  individual  items 
affecting  the  probability  distribution  of  the  indices  (29:8). 
Because  the  indices  of  stochastic  inventories  are  random  vari¬ 
ables,  the  established  objectives  must  be  viewed  in  terms  of 
expected  values  which  are  unlikely  to  be  exactly  met. 

Inventories  with  stochastic  demands  and  leadtimes 
also  require  safety  stock  to  reduce  the  possibility  of  stock¬ 
outs  caused  by  leadtime  demand  forecast  errors.  Safety  stocks 
can  be  defined  as  the  difference  between  the  expected  demand 
for  a  period  and  the  level  of  stocks  held  to  meet  demand  for 
the  period  (31:42).  In  the  continuous  revie\v-  situation,  the 
stocks  are  ordered  when  the  inventory  level  reaches  the  re¬ 
order  point.  The  difference  between  the  reorder  point  and 
the  safety  stock  level  is  the  expected  or  mean  demand  during 
the  leadtime  before  the  next  order  is  received.  The  safety 
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stock,  then,  is  used  to  meet  demands  in  excess  of  the  ex¬ 
pected  demand.  This  excess  demand  represents  the  forecast 
error  of  demand  during  the  leadtime  period.  Because  the  de¬ 
mand  forecast  error  may  exceed  the  level  of  safety  stocks 
held,  stock-outs  may  occur  and  backorders  will  accumulate. 

As  with  the  overall  inventory,  optimum  safety  stock 
levels  are  determined  by  consideration  of  the  costs  and 
benefits  of  stocking  various  levels.  The  same  basic  relation¬ 
ships  exist.  As  the  safety  stock  level  is  increased,  the 
level  of  service  increases,  holding  costs  increase,  and 
ordering  costs  may  decrease  (31:42).  At  some  point,  an  opti¬ 
mum  balance  exists  between  the  desired  stock  level,  service 
level,  and  costs.  In  multiple- item  inventories,  these  effec¬ 
tiveness  measures  can  again  be  represented  by  various 
aggregate  indices  to  indicate  the  performance  over  all  items. 
These  safety  stock  indices  are,  of  course,  based  on  single 
items  with  stochastic  leadtime  demands;  therefore,  the  indices 
are  stochastic  and  can  be  represented  as  probability  distri¬ 
butions  which  are  affected  by  the  leadtime  demand  patterns  of 
the  individual  items.  The  safety  stock  policy  for  a  multiple- 
item  inventory  represents  an  attempt  to  find  an  optimum  bal¬ 
ance  between  the  applicable  costs,  safety  stock  level,  and 
service  level  as  measured  by  stochastic  aggregate  indices. 

The  current  DESC  safety  stock  policy  is  determined  by  the  DLA 
model  which  is  discussed  in  the  next  section. 


Analysis  of  the  DLA  Model 


The  DLA  model  currently  used  by  DESC  to  compute  vari¬ 
able  safety  stock  levels  stems  from  techniques  and  general 
models  developed  by  Hadley  and  Whitin  for  continuous  review 
inventories  with  stochastic  demands.  In  developing  their 
stochastic  models,  Hadley  and  Whitin  assumed  that  "the  process 
generating  the  demands  does  not  change  with  time"  (i.e.,  the 
mean  rate  of  demand  is  constant  over  time)  and  used  the  mini¬ 
mization  of  variable  costs  as  the  criterion  to  determine  the 
optimal  inventory  policy.  They  also  caution  that,  by  defini¬ 
tion,  a  continuous  review  model  assumes  that  an  order  is 
placed  precisely  when  the  reorder  point  is  reached.  This,  in 
turn,  implies  that  the  number  of  units  demanded  per  requisi¬ 
tion  cannot  be  a  random  variable  because  then  it  would  be 
possible  to  overshoot  the  reorder  point  before  an  order  could 
be  placed.  Therefore,  in  cases  where  the  units  per  requisition 
is  a  random  variable,  it  is  inappropriate  to  order  a  fixed 
quantity  each  time  an  order  is  placed.  Hadley  and  Whitin  note 
that  for  practical  purposes  continuous  review  models  can  be 
used  for  inventories  with  a  variable  requisition  size,  pro¬ 
vided  the  probability  of  overshooting  the  reorder  point  is 
very  small.  In  such  cases,  the  optimal  reorder  point  R  and 
order  quantity  Q  can  be  found  with  a  continuous  review  model, 
then  the  actual  order  quantity  is  Q  plus  the  amount  by  which 
the  reorder  point  was  exceeded  (9:Sec  4-1).  This  is  the 
method  used  by  DESC. 


The  DLA  continuous  review  model  is  based  upon  the 
development  by  Hadley  and  Whitin  of  a  model  which  assumes  that 
item  leadtime  demand  is  approximately  normally  distributed. 

They  note  that  for  large  leadtime  demands,  discreteness  can 
be  ignored  and  all  variables  can  be  treated  as  continuous. 

In  addition,  the  continuous  normal  distribution  will  be  a  good 
approximation  for  the  actual  demand  distribution  when  the 
demand  is  sufficiently  large.  The  model  also  assumes  that 
procurement  leadtimes  are  constant,  and  that  the  entire  order 
is  received  at  the  end  of  the  leadtime.  The  model  allows 
backorders  and  uses  the  inventory  position,  rather  than  the 
inventory  on  hand,  to  define  the  reorder  point.  The  objec¬ 
tive  of  the  Hadley  and  Whitin  model  is  to  develop  an  expres¬ 
sion  for  the  average  annual  cost  of  the  inventory  which  should 
be  minimized  to  find  the  optimal  order  quantity  and  reorder 
point.  However,  in  the  course  of  developing  the  cost  func¬ 
tion,  it  was  also  necessary  to  determine  the  expected  number 
of  units  in  a  backorder  status  at  a  random  point  in  time,  given 
the  assumptions  of  the  model  (9:  Sec  4-9j  .  It  is  this  expres¬ 
sion  for  expected  backorders  that  was  used  in  developing  the 
DLA  model  (18:246) . 

In  their  development  of  stochastic  models,  Hadley 
and  Whitin  also  address  the  difficulty  in  determining  the 
cost  of  a  backorder  and  produce  a  model  which  does  not  require 
that  explicit  backorder  costs  be  assigned.  The  objective  of 
this  model  is  to  minimize  the  average  annual  ordering  and 
holding  costs  subject  to  a  constraint  on  the  expected  number 
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of  units  in  a  backorder  status  at  any  point  in  time.  Using 
the  theory  of  Lagrange  multipliers,  Hadley  and  Whitin  show 
that  this  is  equivalent  to  minimizing  (in  simplified  form): 


TVC  =  OC  +  HC  +  AB 


where 


TVC  =  total  variable  cost 
OC  =  ordering  cost 
HC  =  holding  cost 
B  =  expected  number  of  backorders 
\  =  Lagrange  multiplier 

This  general  model  will  determine  the  optimum  order  quantity 
and  reorder  point  for  an  item  and  uniquely  determine  X,  which 
is  the  implicit  cost  of  a  backorder.  Thus,  by  specifying  a 
backorder  constraint,  a  unique  implicit  backorder  cost  is 
determined  (9:  Sec  4-16). 

This  is  precisely  the  same  model,  presented  in  Chap¬ 
ter  1,  which  was  adopted  as  the  DoD  policy  for  determining 
order  quantities  and  safety  stock  levels.  The  specific  order¬ 
ing  and  holding  cost  functions  used  in  the  DoD  model  are  also 
exactly  the  same  functions  developed  in  the  Hadley  and  Whitin 
model  (9:219).  In  addition,  although  the  DoD  model  does  not 
specify  a  particular  demand  distribution,  the  expression 
given  to  determine  requisitions  short  is  simply  a  generalized 
form  of  the  expected  backorders  expression  developed  by  Hadley 
and  Whitin  for  their  model  with  normally  distributed  demands 
(9:193).  The  only  change  made  by  the  DoD  in  adopting  the 
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Hadley  and  Whitin  model  was  to  multiply  the  expression  for 
expected  backorders  by  E/S  (see  page  14) ,  which  has  the  effect 
of  weighting  the  backorders  with  an  essentiality  factor  and 
converting  the  number  of  units  backordered  to  the  number  of 
requisitions  backordered.  Because  the  DoD  model  is  the  same 
model  as  developed  by  Hadley  and  Whitin,  the  same  assumptions 
must  apply,  although  none  are  specified  for  the  model  pre¬ 
sented  in  DoD  Instruction  4140.39.  Both  models  also  share  a 
computational  difficulty  in  that  they  require  evaluation  of 
a  complex  integral  to  determine  the  number  of  expected  back¬ 
orders  . 

The  more  specific  DLA  model  basically  makes  two  ex¬ 
tensions  to  the  DoD  model.  First,  in  developing  the  model 
used  by  DLA,  Presutti  and  Trepp  substituted  for  the  normal 
distribution  a  probability  density  function  which  approximates 
the  normal  distribution,  but  can  be  easily  integrated.  This 
eliminated  the  computational  problem  of  the  DoD  model  by 
allowing  the  expected  number  of  backorders  to  be  explicitly 
defined  as  an  easily  calculated  expression  rather  than  in 
terms  of  a  complex  integral.  Secondly,  Presutti  and  Trepp 
explicitly  define  the  safety  stock  level  in  terms  of  the 
standard  deviation  of  the  leadtime  demand  multiplied  by  some 
safety  factor.  These  two  extensions  to  the  DoD  model  allow 
the  safety  factor  and,  therefore,  the  safety  level  to  be  ex¬ 
plicitly  defined  by  an  easily  calculated  expression.  The 
assumptions  upon  which  Presutti  and  Trepp  based  the  develop¬ 
ment  of  the  DLA  model  are  the  same  as  those  used  by  Hadley 
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and  Whitin.  In  addition,  of  course,  they  assumed  that  the 
substitute  probability  density  function  is  an  adequate  appro¬ 
ximation  for  the  normal  distribution  over  the  range  of 
interest  (18) . 

In  using  the  model  to  calculate  safety  levels,  DLA 
also  assumes  that  the  standard  deviation  of  the  leadtime  de¬ 
mand  equals  1.25  times  the  mean  absolute  deviation  of  the 
leadtime  demand  forecast  errors.  This  assumption  is  valid 
when  the  forecast  errors  are  normally  distributed  (4:242). 

DLA  also  uses  a  constant  holding  cost  rate  for  all  inventory 
items.  The  current  model  used  by  DESC ,  then,  is  based  on 
several  important  assumptions  which  can  be  summarized  as 
follows : 

1.  The  mean  rate  of  demand  for  each  item  remains 
constant  over  time. 

2.  The  number  of  units  per  requisition  is  not  a 
random  variable. 

3.  The  procurement  leadtime  for  each  item  is  constant. 

4.  The  entire  order  is  received  at  the  end  of  the 
leadtime . 

5.  Leadtime  demands  are  large  enough  for  the  demancs 
to  be  approximately  normally  distributed. 

6.  The  substitute  function  adequately  approximates 
the  normal  distribution  over  the  relevant  range. 

7.  The  demand  forecast  errors  are  normally  distributed. 

8.  The  holding  cost  rate  is  constant  for  all  items. 

It  is  difficult  to  determine  the  validity  of  these 


and  other  assumptions  which  may  be  implicit  in  the  model  or 
to  determine  the  aggregate  effect  of  any  deviations.  In  addi¬ 
tion  to  the  question  of  whether  it  is  a  reasonably  accurate 
representation  of  the  DESC  inventory,  the  model  has  some  char¬ 
acteristics  which  should  be  noted.  First,  the  model  allows 
the  customer  service  level  of  the  inventory  to  be  managed  as 
an  aggregate  index  which  is  measured  by  the  total  number  of 
requisitions  short.  Also,  the  order  quantity  in  this  model 
is  completely  independent  of  the  backorder  constraint  which 
allows  the  level  of  customer  service  to  be  adjusted  by  chang¬ 
ing  only  the  safety  stock  levels  and  using  the  same  economic 
order  quantities  (18:249).  In  addition,  for  a  given  service 
level,  or  backorder  constraint,  and  any  set  of  order  quanti¬ 
ties,  this  model  will  minimize  the  cost  of  holding  the  safety 
stock  (18:250).  However,  the  model  may  not  allocate  adequate 
safety  stock  levels  to  high  cost  items.  Since  it  is  less  ex¬ 
pensive  to  procure  lower  cost  items  to  reduce  the  number  of 
backorders  for  the  total  inventory,  high  cost  items  will  have 
lower  safety  levels.  Finally,  use  of  the  model  implicitly 
assumes  that  enough  funds  are  available  to  purchase  the  levels 
of  safety  stock  determined  by  the  backorder  constraint.  When 
funds  are  restricted,  it  is  necessary  to  be  able  to  select  the 
backorder  constraint  on  the  basis  of  the  available  funds. 

Relating  Safety  Stock  Level  to  Service  Level 
In  discussing  restricted  funds,  Hadley  and  Whitin 
"None  of  the  inventory  models  discussed  in  operations 
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research  appears  to  provide  an  adequate  treatment  of  budget 
constraint  [10:152]."  Like  most  models,  the  DLA  model  balances 
the  inventory  operating  costs  with  the  service  level  without 
considering  the  stock  levels  as  represented  by  the  budget  con¬ 
straint.  In  the  DLA  model,  the  service  level,  measured  as 
the  number  of  requisitions  short,  is  determined  solely  by  the 
safety  stock  level,  which  is  measured  in  dollars  invested  in 
safety  stock.  Therefore,  when  the  inventory  managers  have 
limited  funds  to  invest  in  safety  stock  or  want  to  balance 
investments  against  service,  they  must  be  able  to  determine 
the  variable  safety  level  dollars  (VSL$)  required  for  any 
specified  level  of  backorders  (8).  In  establishing  the  model 
to  be  used  to  calculate  safety  levels,  DoD  Instruction  4140.39 
notes  that  the  shortage  parameter  X  can  be  selected  to  both 
control  the  safety  level  and  to  satisfy  other  constraints  such 
as  the  funding  level  (28:2).  This  implies  that  X  might  be 
used  to  relate  VSLS  to  g  and  select  an  appropriate  balance; 
however,  no  attempt  was  made  to  define  this  relationship. 

Several  studies  concerned  with  the  model  (3,  5,  7, 

15)  also  note  that  X  can  be  set  to  control  8  and/or  VSl$. 
However,  again  no  method  is  suggested  about  how  X  could  be 
used  as  a  control  or  how  a  relationship  between  3  and  VSL$ 
might  be  established  using  X.  Deemer  emphasizes  that  X  "is 
only  an  implied  shortage  cost,  the  true  cost  being  unknown 
[5:7]."  Therefore,  a  value  of  X  has  little  meaning  in  itself 
and  is  simply  a  number  which  has  resulted  from  inventory 
policy  decisions.  Since  X  is  clearly  a  function  of  g  (see 


page  14) ,  and  8  is  determined  by  VSL$ ,  X  could  possibly  also 
be  defined  as  a  function  of  VSL$ .  These  two  functions  could 
then  be  equated  to  possibly  derive  an  explicit  expression 
relating  VSL$  to  8.  However,  it  makes  little  sense  to  compli¬ 
cate  the  problem  by  including  an  implied  cost  rather  than 
attempting  to  find  a  direct  relationship  between  VSL$  and  B. 

Presutti  and  Trepp  suggest  that  different  values  of 
X  be  used  to  calculate  the  values  for  S  and  VSL$  from  the  cur¬ 
rent  inventory  data,  then  to  plot  the  values  to  obtain  two 
curves  (18:250).  One  curve  will  relate  X  to  service  level 
and  the  other  to  investment  level.  This  procedure  has  the 
disadvantage,  among  others  to  be  discussed  shortly,  that  it 
also  unnecessarily  relates  8  and  VSL$  to  X  rather  than  directly 
to  each  other.  Another  indirect  approach,  used  by  the  Navy 
on  a  different  form  of  the  DoD  model,  is  to  calculate  service 
levels  based  on  changing  values  of  X  and  iterate  to  a  solution 
which  satisfies  the  funding  constraint  (15:21).  Although  this 
technique  produces  a  solution  which  satisfies  the  funding  con¬ 
straint,  the  constraint  is  necessarily  somewhat  arbitrary 
since  the  service  level  is  not  known  until  after  inventory 
policy  is  calculated.  Both  this  procedure  and  the  one  sug¬ 
gested  by  Presutti  and  Trepp  are  also  undesirable  because  they 
require  that  several  values  of  each  variable  be  calculated. 

This  is  necessary  to  obtain  a  reasonably  accurate  curve  or  to 
iterate  to  a  reasonably  close  solution  and  increases  the 
cost  of  the  analysis. 

DLA  suggests  the  direct  approach  that  only  VSL$  and 


g  be  computed  and  then  plotted  against  each  other  to  derive  a 
curve  as  shown  in  Figure  4  (6:3).  The  basic  form  of  this 
curve  is  generally  known  (4:20;  7:111-17)  and  easily  explained. 
At  zero  VSL$ ,  or  no  safety  stock,  there  is  some  maximum  number 
of  expected  backorders.  As  VSL$  increases,  backorders  decrease 
at  a  diminishing  rate  until  eventually  additional  safety  stock 
has  little  effect  on  decreasing  the  number  of  backorders. 

This  approach  also  requires  that  several  values  of  both  vari¬ 
ables  be  calculated  to  obtain  an  accurate  curve.  In  addition, 
a  problem  shared  by  each  of  these  procedures  is  that  the  rela¬ 
tionships  between  the  variables  are  probably  not  stationary 
and  are  valid  only  for  the  period  in  which  they  are  computed. 
Since  demands  and  prices  are  known  to  change  over  time  (1) , 
and  these  effect  g ,  X ,  and  VSL$  ,  the  relationship  between  6 
and  VSL$  changes  over  time.  Therefore,  to  obtain  the  correct 
curve  requires  that  3  and  VSL$  be  recomputed  from  the  current 
inventory  data  each  time  (at  least  quarterly)  that  the  rela¬ 
tionship  is  required  to  make  a  decision.  Due  to  the  importance 
of  being  able  to  select  3  based  on  available  VSL$ ,  DESC  has 
recently  decided  to  use  the  DLA  suggested  approach  and  calcu¬ 
lates  the  backorders  and  VSL$ ,  based  on  the  total  inventory, 
prior  to  each  quarterly  update. 

The  problems  with  these  procedures  are  that  they  do 
not  provide  the  information  required  by  DESC,  and  presumably 
other  DoD  supply  organizations,  or  they  are  very  costly  be¬ 
cause  they  require  many  calculations  on  a  large  data  base. 

Since  the  information  required  is  the  current  relationship 


Relationship  Between  Expected  Backorders 
and  Safety  Stock  Investment 


between  8  and  VSL$ ,  the  approach  suggested  by  DLA  of  deriving 
the  curve  defining  this  relationship  certainly  provides  the 
information  in  the  most  desirable  form.  Also,  although  this 
curve  is  not  stationary  and  must  be  derived  for  each  period, 
it  is  possible  that  the  basic  form  of  the  curve  is  the  same 
for  all  periods.  This  is  suggested  by  the  rationale  of  dim¬ 
inishing  marginal  returns,  which  explains  the  form  of  the 
curve,  and  the  fact  that  the  same  model  is  used  to  derive  the 
curve  for  each  period.  If  the  basic  shape  of  the  curve  is 
the  same  for  each  period,  it  should  be  possible  to  derive  a 
general  equation  which  defines  the  basic  relationship  between 
3  and  VSL$ .  With  such  an  equation,  fewer  calculations  should 
be  required  to  define  the  curve.  The  calculations  required 
to  find  the  equation  exactly  will,  of  course,  be  determined 
by  the  number  of  unknown  variables  in  the  equation.  However, 
the  appavntly  simple  form  of  the  curve  suggests  that  a  rela¬ 
tively  simple  equation  can  be  devised  which  would  significantly 
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reduce  the  computational  burden  of  finding  the  curve.  One  of 
the  research  objectives  was  to  derive  such  an  equation. 

Even  with  a  general  equation  to  define  the  basic  form 
of  the  curve,  the  present  method  would  require  using  current 
data  from  all  inventory  items  to  find  the  curve  for  the  cur¬ 
rent  period.  To  further  reduce  the  cost  of  finding  the  curve, 
it  is  necessary  to  eliminate  or  reduce  the  amount  of  data  re¬ 
quired.  One  method  to  do  this  is  to  predict  the  relationship 
from  past  aggregate  indices  which  are  related  to  6  and  VSL$. 
DESC  calculates  several  aggregate  indices  each  quarter  to 
measure  various  aspects  of  the  total  inventory.  Since  8  and 
VSL$  are  closely  related  to  some  of  these  indices,  it  should 
be  possible  to  predict  the  relationship  between  B  and  VSL$ , 
if  the  other  indices  can  be  accurately  forecasted.  This  pro¬ 
cedure  would  completely  eliminate  the  need  for  data  on  indivi¬ 
dual  items  and  greatly  reduce  the  number  of  calculations  re¬ 
quired.  Another  method  to  reduce  the  number  of  calculations 
is  to  use  a  sample  of  the  current  data  rather  than  data  for 
all  items.  Thus,  another  research  objective  was  to  determine 
how  accurately  the  relationship  between  B  and  VSLS  could  be 
predicted  using  less  data.  The  method  tested  was  to  predict 
the  relationship  based  on  forecasts  of  the  aggregate  indices. 

Some  Alternative  Models 

Although  the  purpose  of  this  research  was  to  find  a 


possible  solution  in  the  context  of  the  current  DLA  model,  it 
must  be  acknowledged  that  a  better  solution  might  be  obtained 
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by  changing  to  a  different  model.  Recognizing  that  most 
inventory  models  are  "severe  abstractions  of  the  real  world 
[3:2],"  other  models  may  set  safety  levels  as  close  to  the 
optimum  as  the  current  model  and  at  the  same  time  eliminate 
or  reduce  the  problem  of  relating  g  to  VSL$ .  The  current  DLA 
model  is  one  of  a  class  of  literally  hundreds  of  traditional 
cost  models  which,  based  on  various  assumptions,  attempt  to 
minimize  the  sum  of  holding,  ordering,  and  stockout  costs. 
Even  within  this  class,  it  is  questionable  that  the  DLA  model 
is  the  most  appropriate  for  the  situation.  In  fact,  in  dis¬ 
cussing  a  case  where  this  type  of  model  was  applied  to  "a 
military  supply  system  concerned  with  stockage  of  electronic 
components,"  Hadley  and  Whitin  concluded  that  their  model  was 
"entirely  inappropriate.  .  .[in]  a  situation  where  there  was 
a  fixed  annual  procurement  budget  [9:403]."  A  variation  of 
the  conventional  model  which  may  be  more  appropriate  in  this 
situation  is  to  minimize  backorders  subject  to  a  constraint 
on  inventory  investment.  Detailed  models  of  this  form  have 
been  developed  (20,  26)  and  are  a  possible  alternative  to  the 
current  model. 

Another  type  of  model  which  may  be  an  improvement 
over  the  DLA  model  is  some  form  of  dynamic  programming  model. 
The  advantage  of  this  class  of  models  is  that  they  take  into 
consideration  the  changes  which  occur  in  inventory  systems 
over  time.  In  particular,  dynamic  models  are  helpful  when 
demand  distributions,  prices,  and  the  items  carried  in  the 
inventory  are  changing  (9:323).  Since  these  do  change  in  the 


DESC  inventory,  a  dynamic  model  may  provide  more  appropriate 
safety  levels  and  allow  the  relationship  between  B  and  VSL$ 
to  be  determined  more  accurately.  An  important  disadvantage 
of  dynamic  models  is  that  they  are  generally  much  more  complex 
and,  therefore,  more  costly  to  use  than  steady-state  models. 

One  study  which  developed  a  dynamic  programming  model  that 
minimizes  expected  shortages  subject  to  an  investment  constraint 
concluded  that  the  model  was  computationally  impractical  for 
inventories  as  large  as  100,000  items.  However,  the  study 
suggested  that  if  items  could  be  grouped  according  to  demand 
rates  and  costs,  the  model  could  be  used  to  allocate  shortages 
and  funds  to  the  groups  (12:32). 

Another  alternative  is  to  use  a  model  based  upon  mar¬ 
ginal  returns.  As  noted  above,  the  basic  relationship  between 
6  and  VSL$  is  one  of  decreasing  marginal  returns.  As  the 
number  of  backorders  decreases,  the  ratio  of  B/VSL$  increases. 

A  model  based  on  marginal  returns  iteratively  allocates  funds 
to  the  item  which  maximizes  the  desired  performance  measure 
per  dollar  invested  (19:3).  In  the  case  of  backorders, 
safety  level  dollars  are  allocated  iteratively  to  the  item 
with  the  maximum  B/VSL$  ratio.  The  allocation  of  VSL$  con¬ 
tinues  until  the  investment  constraint  is  reached  or  an 
acceptable  backorder  level  is  reached.  A  multi - echelon  model 
of  this  form,  based  on  minimizing  backorders  with  an  invest¬ 
ment  constraint,  was  developed  for  the  Air  Force  base-depot 
supply  system  (21).  Application  of  a  marginal  return  model 
to  base  level  consumable  secondary  items  has  also  been 


studied  (19).  It  appears  that  this  form  of  model  has  not  yet 
been  applied  to  any  problem  as  large  as  the  DESC  inventory, 
perhaps  because  of  the  large  number  of  computations  required 
for  so  many  items.  This  problem,  as  with  the  dynamic  pro¬ 
gramming  model,  might  be  solved  by  grouping  items  and 
allocating  VSL$  among  the  groups. 

One  more  approach  to  the  problem  involves  using  the 
current  model  to  determine  an  initial  inventory  policy  and 
another  model  to  alter  the  policy  if  sufficient  VSL$  are  not 
available.  A  model  has  been  developed  for  the  Army  which 
modifies  inventory  policies  when  insufficient  funds  are  avail¬ 
able  to  buy  the  levels  of  stock  indicated  by  the  current 
policy  (13:4).  This  model  evaluates  the  changes  in  required 
investment  if  a  particular  supply  modification,  such  as  cut¬ 
ting  reorder  points  and  raising  backorders,  is  implemented. 
After  evaluating  several  possible  modifications,  the  policy 
is  chosen  which  meets  the  funding  constraint  and  provides  the 
highest  level  of  service.  While  this  model  was  developed  to 
alleviate  a  problem  essentially  the  same  as  the  one  con¬ 
sidered  here,  it  appears  to  be  an  unnecessarily  indirect  and 
expensive  approach. 


Summary 

The  purpose  of  this  chapter  was  to  thoroughly  review 


the  research  problem  and  to  develop  and  support  the  specific 
research  objectives.  After  discussing  the  basic  theory  of 
managing  safety  stocks  in  mult iple- item  stochastic  inventories, 
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the  theory  underlying  the  DLA  model  was  reviewed.  The  prin¬ 
ciple  assumptions  upon  which  the  model  is  based  were  empha¬ 
sized  and  some  of  the  more  important  strengths  and  weaknesses 
of  the  model  were  noted.  Several  possible  methods  to  relate 
backorders  to  safety  level  dollars  were  discussed  before  pre¬ 
senting  the  approach  selected  for  this  research.  This 
approach  was  to  attempt  two  methods  to  reduce  the  computational 
effort  required  to  derive  a  relationship  between  backorders 
and  safety  level  dollars.  The  first  method  was  to  derive  a 
general  equation  which  defined  the  basic  form  of  the  6  -  VSLS 
curve.  The  second  method  was  to  reduce  the  amount  of  data 
required  by  using  aggregate  indices  from  past  periods  to 
estimate  the  current  relationship  betiveen  S  and  VSL$ . 

Finally,  some  alternative  models  were  presented  which  showed 
that  there  are  other  approaches  to  the  problem  which  could  be 
pursued . 
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CHAPTER  3 


METHODOLOGY 

The  methodology  adopted  for  this  research  is  based 
on  the  assumption  that  a  causal  relationship  exists  between 
VSL$  and  8,  as  explained  in  Chapter  2,  and  that  the  relation¬ 
ship  can  be  predicted  for  future  periods.  In  addition,  it 
was  assumed  that  the  relationship  during  any  period  would  be 
a  function  of  the  state  of  the  total  inventory  as  measured 
by  various  aggregate  indices.  The  rationale  for  this  assump¬ 
tion  is  that  individual  item  variable  safety  levels  (VSL^) 
are  determined  by  both  6  and  various  item  characteristics 
such  as  price,  requisition  size,  and  leadtime  demand  variance. 
Therefore,  the  aggregate  safety  level  measure  VSL$  =  EVSL^C^ 
should  be  a  function  of  both  8  and  some  aggregate  indices  of 
the  item  characteristics.  Based  on  these  assumptions,  the 
overall  approach  was  to  first  develop  models  which  would  de¬ 
fine  VSL$  as  a  function  of  8  and  other  available  indices  for 
any  period.  Then,  to  determine  the  relationship  between  VSL$ 
and  s  for  the  upcoming  period,  the  values  of  the  indices  used 
in  the  models  would  be  forecasted.  The  methods  used  to  develop 
the  models  and  forecast  the  values  of  the  indices  are  explained 
in  this  chapter.  This  is  preceded  by  a  discussion  of  the  data 
used  and  followed  by  the  method  used  to  test  the  models  and 
forecasting  technique. 


Data  Selection  and  Definitions 


Twenty-eight  sets  of  data  were  available  from  DESC 
covering  the  period  since  the  DLA  model  was  implemented. 

Each  set  contains  the  aggregate  indices  produced  from  quar¬ 
terly  updates  which  set  new  safety  levels  for  all  inventory 
items  based  on  the  characteristics  of  the  items  in  the  inven¬ 
tory  at  the  time.  However,  some  quarters  were  represented 
by  two  sets  of  data  because  the  quarterly  update  had  to  be 
rerun.  This  was  necessary  either  because  some  input  data 
was  incorrect  or  because  the  value  of  8  chosen  for  the  ini¬ 
tial  run  produced  an  unacceptable  VSL$  level.  In  these 
cases,  only  the  final  data  set,  which  established  the  safety 
levels  actually  used  for  the  quarter,  were  used  in  the 
analysis.  In  addition,  the  sets  from  the  last  two  periods 
available  were  reserved  to  test  the  predictive  ability  of 
the  models  developed.  The  model  development  and  analysis, 
then,  was  based  on  21  sets  of  historical  data  representing 
the  quarters  from  January  1976  to  December  1980.  This  data 
and  the  data  for  the  two  quarters  used  to  test  the  models 
are  presented  in  Appendix  A.  The  variable  names  used  for 
the  indices  and  their  definitions  are  given  below: 

8  -  The  expected  (average)  number  of  time-weighted 

requisitions  short  at  any  point  in  time 
SC  -  Referred  to  as  the  system  constant,  this  is 
the  sum,  over  all  items,  of  the  unit  price 
times  the  mean  absolute  deviation  of  the  lead- 

time  demand  forecast  error  (IC.MADLT.) 

^  1  1 1 


OLD  SC  -  The  system  constant  from  the  period  immediately 


VSL$ 

DEMANDS 

ITEMS 
0- ITEMS 
FREQ 


preceding  the  period  in  question 
Variable  safety  level  dollars  is  the  sum  of  the 
dollar  value  of  the  item  safety  levels  (IVSL^C^) 
The  value  of  the  forecasted  inventory  annual 
demand  (E4QFD^C^) 

The  number  of  items 

The  number  of  items  with  a  safety  level  of  zero 
The  annual  demand  frequency  is  the  total 
number  of  requisitions  per  year  for  all  items 


Model  Development 

Correlation  and  regression  analysis  were  the  basic 
methods  used  to  develop  and  compare  possible  models.  Correla¬ 
tion  coefficients  were  used  to  measure  the  strength  of  the 
pairwise  relationships  between  VSL$ ,  3,  and  the  other  avail¬ 
able  indices.  The  correlation  statistics  used  were  the 
Pearson  product -moment  correlation  coefficient  r  and  a  signi¬ 
ficance  test  on  r  using  the  Student’s  t  with  N-2  degrees  of 
freedom  for  the  quantity  (17:281): 


This  information  was  then  used  to  help  determine  which  vari¬ 
ables  might  be  added  or  deleted  in  a  model  to  improve  the 
model.  Multiple  regression  analysis  was  used  to  develop 
models  and  was  the  primary  method  used  to  compare  the  models 
The  DLA  model  was  also  examined  and  manipulated  to  determine 
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what  forms  of  regression  models  might  be  suggested  by  the 
structure  of  the  relationships  between  the  variables  in  the 
DLA  model.  Finally,  the  basic  VSL$  -  B  curve  suggested  a 
model  form  which  was  then  further  developed  through  regres¬ 
sion  analysis.  The  model  development  process  eventually  pro¬ 
duced  three  basic  forms  of  models,  which  were  compared  to 
select  the  one  which  most  accurately  defined  the  relationship 
between  VSL$  and  B.  The  models  which  were  developed  are  pre¬ 
sented  below. 

Linear  Multiple  Regressions 

Two  models  were  developed  using  simple  linear  mul¬ 
tiple  regression  techniques.  In  one  model,  VSL$  was  the 
independent  variable,  and  in  the  second  S  was  the  independent 
variable.  These  models  were  developed  using  the  stepwise 
regression  procedure  from  the  Statistical  Package  for  the 
Social  Sciences  (SPSS)  (17).  In  this  procedure,  one  variable 
is  chosen  at  each  step  to  be  added  or  removed  from  the  regres¬ 
sion  equation.  The  selection  process  is  based  on  the  partial 
F- value  of  each  variable  which  indicates  the  amount  each  vari¬ 
able  contributes  to  reducing  the  residual  sum  of  squares.  If 
the  F-value  of  a  variable  falls  below  a  specified  value  during 
the  stepwise  procedure,  it  is  removed  from  the  equation  in 
the  next  step.  If  no  variables  are  to  be  removed  in  a  step, 
the  variable  with  the  highest  F-value  and  not  already  in  the 
equation  is  entered.  This  procedure  continues  until  all  the 
variables  not  in  the  equation  have  F-values  below  the  specified 
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value.  The  F-value  specified  to  enter  and  remove  variables 
from  the  equation  was  2.0  for  all  of  the  SPSS  regressions  per¬ 
formed  in  this  study.  The  form  of  the  two  linear  models  which 
will  be  referred  to  as  Models  1  and  2  was: 

Model  1:  VSL$  =  CQ  +  C1B  +  +  C3X2  +  ...  +  C^X. 

Model  2:  6  =  CQ  +  CjVSLS  +  +  C3X  +  ...  +  CjX. 

where  the  X-  are  the  inventory  indices  and  the  are  the 
regression  coefficients. 

Regressions  Based  on  the 
DLA  ModeT 

The  second  form  of  regression  model  tested  was  derived 
from  the  DLA  variable  safety  level  model.  The  objective  in 
developing  this  form  of  model  was  to  rearrange  and,  where 
possible,  simplify  the  DLA  model  to  obtain  an  expression  which 
represented  a  linear  or  nonlinear  regression  model  with  the 
available  indices  as  the  variables.  The  dependent  variable 
in  the  model  is  VSL$  defined  as: 

VSL$  -  E  VSLiCi 
or 

VSL$  =  E  1.25  kiMADLTiCi 

The  expression  for  k^  (see  page  14)  can  be  rewritten  as: 

k  .  _Llnr. _ Wi _ 

i  a  lZiMADLT.  (l-exp(-/7Qi/l.  25  MADLT^) 

2.56  B  , 

*  Z  C.'MADLT.-1 
i  l 
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or 


k. 

1 


SQ.C. 

ixi  1 


*/2  ln  ^ZiMADLTi(l-exp(-/lQi/1.25  MADLT^) ) 


-  -1=  ln(2.  568)  +  -y=  InCZ  C-MADLT.) 

/2  /2 

This  expression  of  isolates  the  aggregate  indices  from  the 
i^  item  characteristics.  If  X  is  substituted  for  the  term 
in  brackets  containing  all  the  item  characteristics  and  the 
resulting  expression  is  substituted  for  k^  ,  then, 


VSL$  =  1.25  Z  MADLT.Ci(-~)  ln(X) 


♦  1.25  I  MADL^C-e-i)  ln(2.566) 

+  1.25  Z  MADLT.C.  ln(Z  MADLT.C.) 

i  i  /2  1  i 


The  first  term  in  this  expression  is  not  represented  by  any 
available  index,  but  is  a  constant,  Cq ,  in  any  period.  Group¬ 
ing  the  constants  in  the  other  two  terms  and  substituting 
SC  for  Z  MADLT^C^  yields  the  model: 

Model  3:  VSL$  =  Cn  <■  C.SC  ln(B)  ♦  C,SC  ln(SC) 

Because  DESC  uses  the  system  constant  from  the  previous  period 
instead  of  the  last  SC  in  this  expression  to  calculate  safety 
levels,  another  model  of  this  form  was  tested. 


Model  4:  VSL$  =  CQ  +  CjSC  ln(B)  +  C2SC  ln(OLD  SC) 

The  third  model  developed  in  this  group  is: 
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C1(OLD  SC)  In  C  6)  +  C2(OLD  SC)  ln(OLD  SC) 


Model  5:  VSL$  =  CQ  + 

This  model  was  tested  to  determine  how  well  the  old  system 
constant  alone  could  predict  the  VSL$  -  6  relationship. 

Models  3,  4,  and  5  represent  nonlinear  relationships. 
However,  they  are  in  the  form  of  a  linear  regression  model 
when  the  variables  in  each  of  the  second  and  third  terms  are 
combined  and  viewed  as  a  single  regression  variable.  There¬ 
fore,  the  SPSS  linear  stepwise  procedure  was  also  used  to  test 
these  models. 

Model  Based  on  the 
VSL$  -  B  Curve- 

The  development  of  this  model  began  by  plotting  the 
six  VSL$  and  B  values,  provided  by  each  quarterly  update,  for 
several  periods  to  determine  whether  they  actually  fit  the 
curve  hypothesized  in  Chapter  2.  Although  the  values  avail¬ 
able  were  in  a  relatively  small  range,  they  appeared  to  fit 
the  form  of  the  theoretical  curve  as  shown  by  the  examples 
presented  in  Figure  5  .  The  form  of  this  curve  suggested 
that  the  nonlinear  VSL$  -  8  relationship  for  a  single  period 
took  the  form  of  a  power  function  or  logarithmic  function. 

One  form  considered  was 

C2 

VSL$  =  CXB 

which  can  be  transformed  into  the  linear  regression  model 
ln(VSL$ )  =  In (Cl)  +  C2  ln(3) 

Another  model  suggested  by  the  form  of  the  curve  was 
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Figure  5 

Two  Examples  of  Actual  VSL$-B  Curves 
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VSL$  =  Cl  +  C2  ln(B) 


To  help  determine  the  form  of  the  curve,  the  DLA  model  was 
re-examined.  If  Model  3  is  taken  as  a  simplification  of  the 
DLA  model,  it  can  be  further  reduced.  Recognizing  that  the 
SC  is  a  constant  during  any  period.  Model  3  is  transformed 
from 

VSL$  =  CQ  +  C1SC  ln( 2 . 56  0)  +  C2SC  ln(SC) 

to  its  simplest  form 

Model  6:  VSL$  =  CQ  +  Cx  InCS) 

which  is  the  function  suggested  above. 

The  procedure  to  test  Model  6  was  to  first  regress 
the  six  VSL$  values  against  the  ln(B)  values  for  each  period 
to  determine  how  well  the  curve  fit  in  each  period.  These 
regressions  also  produced  the  coefficients  CQ  and  ^  for  each 
period.  To  determine  whether  CQ  and  could  be  determined 
from  the  available  indices,  they  were  regressed  against  these 
variables.  Then  the  regression  equations  for  CQ  and  ^  were 
substituted  into  Model  6  and  the  resulting  regression  model 
was  tested  using  the  SPSS  stepwise  procedure  which  provides 
statistics  indicating  the  accuracy  of  a  regression  model. 
These  were  also  used  to  compare  all  of  the  models  developed. 

Forecasting  Methods 

Although  the  regression  models  indicate  the  relation¬ 
ships  between  VSL$ ,  8,  and  the  other  indices  in  any  period, 
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the  values  of  the  indices  are  not  known  for  the  upcoming 
quarter  which  is  the  period  of  interest.  Therefore,  some 
method  must  be  used  to  determine  the  new  values  of  the  indices 
which  are  included  in  the  regression  models  before  the  quar¬ 
terly  update  is  run.  Then  these  values  can  be  included  in 
the  best  regression  model  to  determine  the  new  VSL$  -  6  curve 
from  which  the  appropriate  8  can  be  selected  to  set  safety 
levels  in  the  quarterly  update.  The  method  selected  was  to 
test  some  basic  forecasting  techniques  on  each  index  and  to 
select  the  most  accurate  technique  for  each. 


Evaluating  Forecasting  Methods 


The  accuracy  of  a  forecasting  technique  can  be  ex¬ 
pressed  in  various  terms  which  are  measures  of  the  average 
amount  of  error  the  technique  produces  for  n  forecasts,  where 


the  error  in  period  t  is  e  »  actual  value  -  forecast  value. 
One  term  is  the  mean  squared  error, 


MSE 


Squaring  the  error  terms  prevents  positive  and  negative  errors 
from  cancelling  each  other  and,  thus,  results  in  a  measure 
of  the  magnitude  of  the  error  expected  in  any  period.  However, 
the  MSE  is  measured  in  the  same  units  as  the  index  being 

forecast.  Therefore,  when  the  index  value  varies  signifi- 

2 

cantly  over  time,  the  e  in  a  period  with  a  high  value  will 

2 

be  greater  than  the  e  ^  in  a  period  with  the  same  relative 
error,  but  a  low  index  value.  The  MSE  indicates  the 
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average  error  over  the  periods  measured,  but  the  error  in  a 
particular  period  is  biased  by  the  index  level.  Another 
measure,  expressed  in  terms  of  the  actual  error  rather  than 
the  squared  error,  is  the  mean  absolute  deviation: 

1  n 

MAD  =  -  Z  I  e,  I 
n  t  =  l  1 

Like  the  MSE ,  the  MAD  is  measured  in  the  units  of  the  index 
being  forecast  and,  therefore,  will  be  biased  by  the  index 
level . 

To  avoid  this  problem,  some  measures  of  forecasting 
accuracy  express  the  error  as  a  percentage  of  the  index  value. 
One  such  measure  is  the  mean  percent  error: 


MPE 


100et 


where  is  the  actual  value  in  period  t.  The  MPE  allows 
positive  and  negative  errors  to  cancel  each  other.  While 
this  may  be  useful  in  measuring  the  cumulative  effect  of  the 
errors,  it  does  not  measure  the  expected  error  for  a  specific 
period.  Another  measure  which  expresses  the  error  as  a  per¬ 
centage  is  the  mean  absolute  percent  error: 


MAPE 


,  n  100 | e .  j 

i  £  - — 

n  t=i  Xt 


This  is  a  measure  of  the  expected  amount  of  error  produced 
by  a  forecasting  technique  which  is  not  biased  by  the  level 
of  the  index  being  forecast.  Therefore,  the  MAPE  was  selected 


as  the  most  appropriate  measure  to  compare  the  accuracy  of 
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the  forecasting  techniques  which  were  tested. 


Moving  Averages 

Two  of  the  techniques  tested  on  each  index  in  the 
regression  models  were  simple  moving  averages  and  double 
moving  averages.  The  method  of  simple  moving  averages  fore¬ 
casts  the  value  of  a  variable  by  using  the  average  value 
over  past  periods  for  the  upcoming  period.  The  numbeT  of 
periods  averaged  can  range  from  one,  which  produces  the  so- 
called  naive  forecast,  to  any  higher  number,  with  the 
smoothing  effect  increasing  as  the  number  of  periods  increases. 
The  smoothing  effect  reduces  the  effect  on  the  forecast  value 
of  random  variations  from  one  period  to  the  next.  The  tech¬ 
nique  of  simple  moving  averages  can  be  represented  as  follows: 


where 


S 


t  +  1 


X  + 
At- 1 


+ 


-n  +  1 


) 


St  =  forecast  for  period  t 
X  =  actual  value  for  period  t 
n  =  number  of  periods  included  in  the  average 


Simple  moving  averages  can  produce  useful  short-range  fore¬ 
casts  when  the  variable  tends  to  vary  about  a  fixed  value. 
However,  the  technique  is  not  generally  suitable  when  the 
data  exhibits  trend,  seasonal,  or  cyclical  patterns  (30:34-35). 

The  technique  of  double  moving  averages  is  generally 
more  accurate  when  these  more  complicated  patterns  exist. 
Essentially,  this  method  recognises  that  the  simple  moving 


n  [■Ji'llilltWl  fell 
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average  will  lag  behind  trended  data  by  some  amount.  In 
addition,  a  moving  average  of  the  simple  moving  average  will 
lag  behind  the  simple  moving  average  by  about  the  same  amount. 
Thus,  by  adding  the  difference  between  the  simple  moving 
average  and  the  double  moving  average  back  to  the  simple 
moving  average,  a  more  accurate  forecast  is  produced.  An 
n-period  double  moving  average  can  be  summarized  as  follows: 


S" 


t  +  1 


+  n  +  1 
Vl  n-1 


(st+: 


*  s’t+i) 


S'  +,  is  the  moving  average  of  the  simple  moving  average 
St+1  and  S”  ^  is  the  double  moving  average  forecast.  The 
term  n+l/n-1  is  an  additional  adjustment  which  has  been  found 
to  make  the  technique  more  accurate.  An  important  limitation 
of  this  technique  is  that  2N  periods  of  data  must  be  used,  or 
twice  the  data  required  for  simple  averages  (30:41-44).  Both 
moving  average  techniques  also  require  that  the  best  number 
of  periods  to  average  be  determined.  The  method  used  here 
was  to  calculate  the  MAPE  using  several  different  numbers  of 
periods.  Then  the  number  selected  was  that  which  gave  the 
lowest  MAPE. 


Exponential  Smoothing 

Two  other  techniques  tested  on  the  indices  were  single 
and  double  exponential  smoothing.  An  advantage  that  these 
methods  have  over  moving  averages  is  that  data  from  different 
periods  can  be  given  different  weights.  Single  exponential 
smoothing  gives  some  selected  weight  a  to  the  most  recent 
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observed  value  and  exponent lal ly  decreasing  weights  to  older 
values.  The  general  relationship  is: 

St+1  =  aXt  +  a(l-a)Xtl  +  a(l-a)2Xt  ^  +  a(l-a)5Xt5  +  ... 
This  expression  can  be  rewritten  as: 


S 


t+  1 


aXt  +  ( 1 -  a ) S t 


While  single  exponential  smoothing  can  provide  slightly  better 
forecasts  than  simple  moving  averages  by  given  greater  weight 
to  the  most  recent  data,  it  is  also  generally  not  suitable 
for  data  which  has  a  pattern  of  basic  changes  over  time  (30: 
36-39).  However,  as  with  moving  averages,  double  smoothing 
can  provide  more  accurate  forecasts  for  data  which  has  a 
trend  pattern.  The  same  concept  underlying  the  double  moving 
average  explains  the  double  exponential  smoothing  technique, 
which  can  be  represented  as: 


S" 


t+1 


=  S 


t*l 


1  -a  *■  St  +  1 


SW 


S"t+^  is  the  double  exponential  smoothing  forecast  and  S't+^ 
is  the  single  exponentially  smoothed  Again,  is  an 

additional  adjustment  which  improves  the  accuracy  of  the 
forecast  (30:44-47).  Both  exponential  smoothing  methods  re¬ 
quire  that  the  best  value  of  a  be  determined  and.  this  was 
again  accomplished  by  calculating  the  MAPE  for  several 
values.  Thus,  a  was  selected  based  on  the  lowest  MAPE. 
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Testing  the  Models 

As  mentioned  earlier,  the  two  most  recent  sets  of 
data  were  reserved  to  test  the  models.  After  the  models 
were  developed  and  compared  using  statistical  measures,  the 
index  values  were  forecast  for  these  two  periods  using  the 
best  forecasting  technique  for  each  index.  Then  the  forecast 
values  were  included  in  the  models  to  calculate  VSL$  based 
on  the  same  values  of  6  actually  used  in  the  two  periods. 

The  difference  between  the  actual  VSL$  in  the  period  and  the 
value  of  VSL$  calculated  from  a  model  represents  the  predic¬ 
tion  error  that  would  have  resulted  using  the  model  in  that 
period.  These  results  were  used  as  a  final  comparison  of  how 
well  the  models  may  be  able  to  predict  the  VSL$  -  S  relation¬ 
ship  for  future  periods. 


Summary 

This  chapter  defined  the  variables  in  the  21  sets  of 
data  used  to  develop  the  regression  models.  The  six  models 
were  developed  by  examining  the  theoretical  relationships 
between  the  variables  in  the  DLA  model  and  using  the  SPSS 
stepwise  regression  procedure  to  select  the  most  significant 
variables.  The  SPSS  statistics  were  also  used  as  one  means 
of  evaluating  how  well  the  models  might  predict  the  VSLS  -  8 
relationships.  Four  forecasting  methods  were  then  teL*ed  on 
each  index  required  in  a  regression  model  and  the  best 
technique  for  a  specific  index  was  selected  based  on  the 
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lowest  MAPE.  Finally,  the  models  and  forecasting  techniques 
were  combined  to  predict  a  VSL$  value  for  the  most  recent  two 
periods  and  the  forecast  errors  were  used  to  evaluate  the 
predictive  value  of  the  models. 
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CHAPTER  4 


RESULTS  AND  ANALYSIS 

This  chapter  presents  the  results  obtained  by  attempt¬ 
ing  to  solve  the  research  problem  using  the  methodology  pre¬ 
sented  in  Chapter  3.  First,  the  results  of  the  regression 
analysis  on  the  six  VSL$  -  8  models  are  presented  along  with 
a  discussion  of  the  strengths  and  weaknesses  of  each  model 
based  on  the  regression  statistics.  Then,  the  results  from 
applying  each  forecasting  technique  to  each  of  the  variables 
in  the  regression  equations  is  presented  and  the  best  of  the 
tested  methods  is  selected  for  each  variable.  Finally,  the 
predictive  accuracy  of  each  model  is  demonstrated  by  using 
the  models  to  forecast  the  VSL$  -  S  relationship  for  two 
periods . 


Regression  Analysis 

The  results  of  the  regression  analysis  for  each  group 
of  models  is  presented  first.  This  is  followed  by  compari¬ 
son  of  the  results  for  the  different  models. 

Linear  Models 

The  potential  independent  variables  for  these  two 
models  were  the  available  inventory  indices.  The  correlation 
coefficients  (r)  and  their  significance  (P)  are  given  in 
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TABLE  1 


Correlation  Coefficients  for  the  Linear  Models 


8 

SC 

OLD  SC 

VSLS 

DEMANDS 

ITEMS  0- ITEMS 

FREQ 

(r) 

(P) 

.0791 
.  367 

.4534 

.019 

.  3000 
.093 

.  3833 
.043 

.  5607 
.004 

.9791  .2002 

.001  .192 

0- ITEMS 
(r) 

(P) 

.6205 

.001 

.3832 

.043 

.4308 

.026 

.  2066 
.184 

.  2882 
.103 

.2679 

.120 

ITEMS 

(r) 

(P) 

.1340 
.  281 

.4850 

.013 

.3376 

.067 

.  4063 
.034 

.  5671 
.004 

DEMANDS 

(r) 

(P) 

.  1999 
.  192 

.9753 

.001 

.8146 

.001 

.9574 
.  001 

SL 

(r) 

(P) 

-  .2248 
.  164 

-  .0273 
.453 

-  .0238 
.459 

.1666 
.  235 

VSLS 

(r) 

(P) 

.1343 
.  281 

.9587 

.001 

.7955 

.001 

OLD  SC 
(r) 

(P) 

.6217 

.001 

.8697 

.001 

SC 

(r) 

(P) 

.  3270 
.074 

Table  1.  This  shows  that  S  does  not  have  a  strong  relation¬ 
ship  with  any  individual  variable.  The  highest  correlation 
is  with  the  OLD  SC,  which  explains  only  about  39  percent 
(the  value  of  r  )  of  the  variation  in  B.  VSL$ ,  however,  is 
strongly  related  with  both  SC  and  DEMANDS.  The  correlation 
between  VSLS  and  OLD  SC  is  also  fairly  strong  and  the 
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significance  of  the  three  correlations  is  high.  Therefore, 
a  fairly  accurate  linear  model  was  expected  using  VSL$  as 
the  dependent  variable,  while  Model  2,  with  8  as  the  depend¬ 
ent  variable,  was  not  expected  to  fit  the  data  well. 

The  regression  equation  derived  for  Model  1  was: 

VSL$  =  9.2817  +  .  589 5SC  -  .00048  +  .2084  (OLD  SC) 

Table  2  gives  a  summary  of  the  regression  analysis  for 
Model  1 . 


TABLE  2 

Regression  Summary  for  Model  1 


Step 

Variable 

Entered 

r2 

Overall  F 

Significance 

1 

SC 

.91902 

215.619 

<  .001 

2 

8 

.95499 

190.956 

<  .001 

3 

OLD  SC 

.96244 

145.195 

<  .001 

At  step  3,  the  F  values  of  each  of  the  variables  not  already 
in  the  equation  fell  below  2.0  and,  therefore,  they  were  not 
entered.  DEMANDS  may  have  been  expected  to  enter  the  equa¬ 
tion  because  of  its  high  correlation  with  VSL$ .  ■However,  SC 
entered  the  equation  first  and  it  is  also  highly  correlated 
with  DEMANDS.  Therefore,  DEMANDS  could  explain  little  of 
the  remaining  variation  in  VSL$  which  caused  the  partial  F 
value  for  DEMANDS  to  be  lower  than  the  partial  F  for  both  8 
and  OLD  SC.  An  analysis  of  the  coefficients  in  this  model 
is  shown  in  Table  3. 
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TABLE  3 

Analysis  of  Model  1  Coefficients 


Coefficient 

Standard 

Error 

Partial  F 

Significance 

co 

2.997 

9.591 

.  007 

C1 

0.087 

53.413 

<  .  001 

C2 

0.00009 

16.929 

.  001 

C3 

0.114 

3.371 

.084 

Table  3  shows  that  the  coefficients  have  a  high  significance 

? 

except  for  C^,  the  coefficient  for  OLD  SC.  The  high  r  and 
high  overall  significance  of  Model  1  indicate  that  the  model 
fits  the  data  quite  well.  However,  the  low  significance  of 
C,  reduces  confidence  in  this  model  somewhat  because  this 
indicates  that  old  SC  may  not  have  an  important  influence  on 
VSLS  after  accounting  for  the  influence  of  g  and  SC. 

The  regression  analysis  for  Model  2  produced  the 
following  equation: 

6  =  -4888. 52S6  +  710.296  (OLD  SC) 

-  632.226VSLS  +  1.364(0  ITEMS) 

Table  4  provides  a  summary  of  the  regression  analysis  for 
this  model.  As  expected,  by  the  low  correlation  coefficients, 
Model  2  does  not  fit  the  data  very  well.  Although  the  over- 
all  significance  of  the  model  is  high,  the  r  is  relatively 
low.  This  indicates  a  high  probability  that  the  independent 
variables  contribute  to  determining  the  value  of  VSL$ ,  but 
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TABLE  4 


Regression  Summary  for  Model  2 


TABLE  5 

Analysis  of  Model  2  Coefficients 


Coefficient 

Standard 

Error 

Partial  F 

Significance 

co 

9068.121 

0.290 

.  597 

C1 

118  .  289 

36.057 

<  .001 

C2 

132.546 

22.752 

<  .001 

C3 

0.586 

5.424 

.032 

Step 

Var iab le 
Entered 

r2 

Overall  F 

Significance 

1 

OLD  SC 

. 38644 

11.972 

.003 

2 

VSL$ 

. 74008 

25.626 

<  .001 

3 

0-ITEMS 

. 80295 

23.091 

<  .001 

that  only  about  80  percent  of  the  VSL$  variation  is  explained 
by  the  model.  Table  5  also  shows  that  the  significance  of 
Cj ,  the  coefficient  of  0-ITEMS,  is  not  especially  high,  which 
also  reduces  confidence  in  the  model.  Overall,  Model  2  showed 
little  promise  of  defining  the  VSL$  -  8  relationship  with  the 
desired  accuracy  and  was  rejected  after  examining  the  results 
of  the  regression  analysis. 

Regressions  Based  on  the 
DLA  Model 


The  first  regression  model  derived  from  analyzing  the 
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DLA  safety  level  model  was  Model  3: 

VSL$  =  CQ  +  C:SC  ln(6)  ♦  C2SC  ln(SC) 

This  model  uses  the  current  system  constant  in  all  terms. 

The  correlation  coefficients  for  the  terms  in  this  model  are 
shown  in  Table  6  along  with  the  coefficients  which  apply  to 
Models  4  and  5.  This  table  only  gives  the  correlations  bet¬ 
ween  variables  which  are  in  the  same  model,  as  the  other  rela 
tionships  are  of  no  interest. 

TABLE  6 


Correlation  Coefficients  for  Models  3,  4,  and 

5 

_ 

VSL$  SC  In  (OLD  SC)  ^^(g^  SC 

In ( SC) 

SC  ln(g) 

(r) 

.9365  .9969 

9948 

(P) 

.001  .001 

001 

SC  ln(OLD  SC) 

(r) 

.9498 

(P) 

.001 

(OLD  SC) In ( 6) 

(r) 

.  7  562 

(P) 

.001 

(OLD  SC) In (OLD 

SC) 

(r) 

.8032  .9960 

(P) 

.001  .001 

SC  In (SC) 

(r) 

.9612 

(P) 

.001 

Because  the  term  SC  ln(SC)  has  a  higher  correlation 
with  VSL$ ,  it  would  normally  be  the  first  variable  entered 


in  the  regression  equation  for  the  model  by  the  SPSS  stepwise 
subprogram.  However,  this  may  have  prevented  SC  ln(0)  from 
entering  the  equation  because  of  its  high  correlation  with 
SC  ln(SC)  .  Since  0  is  required  in  the  equation  to  define 
the  VSL$  -  0  relationship,  SC  ln(0)  was  forced  into  the 
regression  model  first.  Surprisingly,  even  with  the  high 
correlation  between  SC  ln(0)  and  SC  ln(SC) ,  the  partial  F 
value  for  SC  ln(SC)  after  SC  ln(0)  entered  the  equation  was 
39.43.  Therefore,  SC  ln(SC)  was  also  entered  and  the  result¬ 
ing  regression  model  was 

VSL$  =  27.067  -  . 1 24SC  ln(0)  +  .382SC  ln(SC) 

A  summary  of  the  regression  analysis  for  Model  3  is  given  in 
Table  7  and  an  analysis  of  the  model  coefficients  is  shown  in 
Table  8. 

TABLE  7 


Regression  Summary  for  Model  3 


Step 

Variable 

Entered 

r2 

Overall  F 

Significance 

1 

SC 

ln(3) 

.87712 

135. 616 

<  .001 

2 

SC 

In (SC) 

.96148 

224.669 

<  .001 

2 

This  model  fits  the  data  very  well  as  shown  by  the  high  r  . 

The  high  overall  F  value  shows  that  the  probability  that 
VSL$  is  not  related  to  the  variables  in  the  model  is  extremely 
low.  Likewise,  the  analysis  of  the  coefficients  shows  that 
they  are  highly  significant  with  relatively  small  standard 
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TABLE  8 


Analysis  of  Model  3  Coefficients 


Coeff ic ient 

Standard 

Error 

Partial  F 

S igni f icance 

co 

0. 0296 

17 . S49 

.  001 

C1 

0.0609 

39.429 

<  .001 

C2 

4.0967 

43.652 

<  .  001 

deviations.  These  results  indicate  that  Model  3  may  provide 
good  predictions  of  VSL$  with  less  than  four  percent  of  the 
variation  in  VSL$  unexplained  by  the  model. 

Model  4  is  the  same  as  Model  3,  except  that  ln(OLD  SC) 
is  substituted  for  ln(SC) .  The  form  of  the  model  is: 

VSL$  =  CQ  +  CXSC  ln(8)  +  C7SC  ln(OLD  SC) 

SC  ln(8)  was  again  forced  into  the  regression  equation  first 
to  ensure  that  3  was  in  the  model.  The  regression  equation 
resulting  from  Model  4  was: 

VSL$  =  21.410  -  0.108SC  ln(B)  +  '  '63SC  ln(OLD  SC) 

The  regression  summary  for  this  cdel  u  shown  in  Table  9, 

and  the  analysis  of  the  coefficients  is  given  in  Table  10. 

Although  the  significance  of  this  model  is  very  high,  the 
2 

value  of  r  leaves  over  eight  percent  of  the  variation  in 
VSL$  unexplained.  In  addition,  the  standard  deviations  of 
the  coefficients  are  relatively  large  and  the  significance 
of  C^  is  not  high.  This  model  is  inferior  to  Model  3  in 
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TABLE  9 


Regression  Summary  for  Model  4 


Step 

Variable 

Entered 

r2 

Overall  F 

Significance 

1 

SC  ln(8) 

.87712 

135.616 

<  .  001 

2 

SC  ln(OLD 

SC) .91901 

.02.128 

<  .  001 

TABLE 

10 

Analysis  of  Model 

4  Coefficients 

Standard 

Loerricient  „ 

Error 

Partial  F 

Significance 

co 

6. 0745 

12.422 

O 

o 

C1 

0 . 0555 

3.7753 

.068 

C2 

0.1189 

9.3119 

.  007 

these  respects  and  still  requires  that  the  current  system 
constant  be  known.  Therefore,  this  model  was  rejected  at 
this  point  in  favor  of  Model  3. 

The  purpose  of  Model  5  was  to  substitute  OLD  SC  for 
SC  to  avoid  the  requirement  to  know  the  current  system  con¬ 
stant.  This  model  had  the  general  form 

VSL$  =  CQ  +  Cj.f'OLD  SC)  ln(8)  +  C2(OLD  SC)  ln(OLD  SC) 

Once  again,  it  was  necessary  to  force  the  first  variable  in 
the  equation  first  because  it  had  a  lower  correlation  with 
VSL$  than  the  second  variable.  The  regression  equation  for 
this  model  was: 


64 


Regression  Summary  tor  Model  5 


Step 

Variable 

Entered 

2 

r 

Overall  F 

Signi f icance 

1 

2 

(OLD  SC) In ( 6 ) 

(OLD  SC) In ( OLD  SC) 

.57178 

.88400 

25.369 

68,  585 

<  .001 

TABLE  12 

Analysis  of  Model  5  Coefficients 


Coefficient 

Standard 

Error 

Partial  F 

Signi f icance 

C0 

7 . 2929 

48  ,  306 

<  .001 

C1 

0.0650 

37.072 

<  .001 

C  2 

0.1419 

48.448 

<  .  001 

VSLS  =  50.687  -  0.402(OLD  SC)ln(S) 

+  0.988 (OLD  SC) In (OLD  SC) 

The  regression  summary  and  analysis  of  the  coefficients  are 
given  in  Tables  11  and  12.  This  model  has  a  high  overall 
significance  and  there  is  no  problem  with  the  significance 
or  standard  deviation  of  any  coefficient.  However,  like 
Model  4,  the  value  of  r  is  low,  meaning  that  a  large  per- 
centage  of  the  variation  in  VSL$  is  left  unexplained  by  the 
;  model.  This  appears  to  indicate  that  the  current  system 

constant  is  required  to  accurately  relate  VSL$  and  8. 
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Model  Based  on  the 
VSLS  -  3  Curve 

The  first  step  in  analysing  Model  6  was  to  test  the 
relationship  between  VSL$  and  ln(3)  for  each  period  using 
the  regression  model: 

VSL$  =  C0  +  Cl  In ( S) 

These  results  are  given  in  Appendix  B,  which  shows  the  value 
of  C  ,  C^,  and  the  r"  for  each  period.  The  average  r  for 
the  21  periods  is  greater  than  0.998,  showing  that  the  data 
fits  this  model  almost  perfectly.  The  next  step  in  the 
analysis  of  this  model  was  to  determine  whether  and 
were  dependent  on  the  available  indices.  The  correlations 
between  CQ,  C^,  and  the  indices  are  given  in  Table  13.  The 
correlations  between  the  other  variables  were  included  in 
Table  1. 


TABLE  13 

Correlation  Coefficients  for  Model  6 


SC 

OLD  SC 

DEMANDS 

co 

(r) 

.  9673 

.9551 

.9152 

(P) 

.001 

.001 

.  001 

C1 

(r) 

.9749 

.9085 

.9362 

(P) 

.  001 

.001 

.001 

ITEMS 

0- ITEMS 

FREQ 

C1 

co 

(r) 

.3547 

.4431 

.3131 

.  9939 

(P) 

.057 

.022 

.084 

.  001 

C1 

(r) 

.  3559 

.4076 

.  3216 

(P) 

.057 

.033 

.078 
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Cq  and  C ^  were  then  regressed  against  the  indices  to  derive 
the  following  equations: 

CQ  =  51.399  +  1.988SC  «■  1.086(OLD  SC)  -  2 1  .  7 97 FREQ 

Cx  •  9.125  +  0.451SC  +  0.121(OLD  SC)  -  4  539 FREQ 

The  regression  summaries  for  these  two  equations  are  given 
in  Tables  14  and  15,  which  show  that  Cq  and  C^  can  be  deter¬ 
mined  very  accurately  given  the  necessary  indices.  The 
analysis  of  the  coefficients  is  not  presented  because  the 
purpose  was  not  to  be  able  to  actually  determine  Cq  and  C^ , 
but  to  find  the  indices  which  should  be  substituted  for  them 
in  Model  6.  The  coefficients  in  the  final  form  of  Model  6 
were  derived  by  a  separate  regression  analysis. 


TABLE  14 

Regression  Summary  for  Cq 


Step 

Variable 

Entered 

r  " 

Overall  F 

Significance 

1 

SC 

.93560 

276.046 

<  .001 

U 

OLD  SC 

.97175 

309.634 

<  .  001 

3 

FREQ 

.98223 

313. 271 

<  .  001 

TABLE 

15 

Regression  Summary  for  C^ 

Step 

Variable 

Entered 

7 

r  ~ 

Overall  F 

Signif icance 

1 

SC 

. 95040 

364.091 

<  .001 

7 

FREQ 

.96866 

2. 7  8  . 1 9  5 

<  .001 

3 

OLD  SC 

.97781 

249.697 

<  .  001 
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Substituting  the  appropriate  indices  for  Cq  and 
in  Model  6  gives  the  following  general  model  : 

VSL$  =  (SC+(OLD  SC) +FREQ)  +  (SC+(OLD  SC) +FREQ) In (B) 

This  was  transformed  into  the  following  regression  model: 

VSL$  =  CQ  +  C1SC  +  C ^ (OLD  SC)  +  C^FREQ  +  C4SC  In  (8) 

+  Cs(OLD  SC)  In ( 3)  +  C6FREQ  ln(S) 

This  was  the  final  form  for  Model  6.  The  SPSS  regression 
analysis  produced  the  following  equation: 

VSL$  =  2S.S81  +  7.27 ISC  -  103.019FREQ  -  .621SC  ln(3) 

+  9.106FREQ  In (8) 

The  regression  summary  and  the  analysis  of  the  regression 
coefficients  for  Model  6  are  presented  in  Tables  16  and  17. 

TABLE  16 


Regression  Summary  for  Model  6 


Step 

Variab le 
Entered 

2 

r 

Overall  F 

Significance 

1 

SC 

.91902 

215.619 

< 

.001 

- 

SC  ln(6) 

.95322 

1  83.385 

< 

.001 

- 

FREQ 

.96300 

147.506 

< 

.001 

4 

FREQ  In  (3) 

.97400 

149.822 

< 

.001 

7 

The  high  r  shows  that  this  model  fits  the  data  very  well, 
with  less  than  three  percent  of  the  VSL$  variation  unexplained. 
Model  6  also  has  a  very  high  overall  significance  and  the 
significance  of  the  individual  coefficients  is  relatively 
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TABLE  17 


Analysis  of  Model  6  Coefficients 


Coefficient 

Standard 

Error 

Partial  F 

Significance 

co 

7.9  708 

10.300 

.  005 

C1 

1.9639 

13.706 

.002 

C-, 

:> 

37.5957 

7  .  509 

.  015 

C4 

0.1870 

11.031 

.004 

C  6 

3.5017 

6.763 

.  019 

high.  Overall 

,  this  model 

could  be  expected 

to  determine 

VSL$  very  accurately,  given  the  values  of  the  independent 
variables . 

Summary  of  Regression  Analyses 

Model  1  was  the  better  of  the  two  linear  models,  with 

2 

an  r  of  better  than  .962.  Model  2  was  rejected  because  it 
fit  the  data  much  worse  than  Model  1.  Model  3  was  the  best 

7 

of  those  derived  from  the  DLA  model,  with  an  r  higher  than 
.961.  Both  Models  4  and  5  were  rejected  because  of  their 
comparatively  low  r  values.  Model  6  had  the  highest  r  of 
.974.  Of  the  three  models  which  had  a  close  fit  with  the 
data,  Model  3  had  the  most  significant  coefficients.  The 
significance  of  one  coefficient  in  Model  1  was  relatively 
low  and  all  the  coefficients  in  Model  6  were  less  significant 
than  those  in  Model  3,  although  the  significance  of  each  was 
still  fairly  high.  Each  of  these  three  were  considered  good 
models,  which  were  further  tested  by  using  forecasts  of  the 
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variables  to  predict  VSL$ . 


Forecasting  Results 

Models  1,  3,  and  6  have  a  total  of  four  independent 
variables  amongst  them:  SC,  g,  OLD  SC,  and  FREQ.  However, 

OLD  SC  is  always  known  and  g  will  be  set  based  on  its  rela¬ 
tionship  with  VSL$ .  Therefore,  it  was  necessary  to  forecast 
only  the  two  unknown  variables,  SC  and  FREQ.  The  initial 
forecasting  analysis  included  data  from  all  of  the  available 
21  periods. 

Figure  6  shows  the  MAPEs  which  result  from  using 
single  exponential  smoothing  with  varying  values  of  a  to  fore¬ 
cast  SC.  Figure  7  shows  the  values  of  MAPE  using  double  ex¬ 
ponential  smoothing  to  forecast  SC.  Figure  8  gives  the  MAPE 
values  which  resulted  from  forecasting  SC  using  single  moving 
averages  and  Figure  9  gives  the  MAPEs  resulting  from  using 
double  moving  averages.  In  these  figures,  the  best  MAPE  value 
for  each  method  is  marked  by  an  asterisk,  while  the  best  over¬ 
all  MAPE  for  the  system  constant  is  marked  with  a  double 
asterisk.  Figures  10  through  13  give  the  same  results  for 
FREQ.  These  figures  show  that  the  lowest  MAPE  for  the  system 
constant  resulted  from  using  double  exponential  smoothing  and 
an  of  0.60.  The  lowest  MAPE  for  FREQ  came  from  using  a 
double  moving  average  over  ten  periods. 

The  forecasts  which  result  from  using  these  two 
forecasting  models  are  plotted  against  the  actual  values  in 
Figures  14  and  15.  The  period  numbers  correspond  with  those 
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in  Appendix  A.  Figure  14  shows  that  the  forecasts  for  SC, 
using  double  exponential  smoothing  with  a  =  .6,  follow  the 
actual  values  very  closely  except  during  periods  six  and 
eight.  The  large  errors  in  these  two  periods  were  obviously 
caused  by  the  extreme  drop  in  periods  six  and  seven,  followed 
by  the  sharp  increase  in  period  eight.  According  to  DESC,  the 
very  low  SC  in  periods  six  and  seven  were  the  result  of  a  pol¬ 
icy  decision  which  was  expected  to  significantly  lower  the  SC 
for  these  periods  (1).  Because  of  this,  these  two  periods 
were  regarded  as  distortions  of  the  normal  SC  trend,  which 
could  influence  the  selection  of  the  best  forecasting  method 
or  the  best  a  or  number  of  periods.  In  Figure  15,  FREQ  shows 
the  same  distortion  in  these  two  periods,  which  again  was 
caused  by  the  policy  change.  Because  the  two  periods  did  not 
represent  the  normal  operating  policy  and  they  significantly 
differ  from  the  other  periods,  it  was  decided  that  they  should 
not  be  used  in  comparing  the  forecasting  methods. 

Therefore,  each  of  the  methods  was  again  tested  on 
both  SC  and  FREQ,  using  only  the  data  from  periods  8  through 
21.  Figures  16  through  19  show  the  new  MAPEs  for  SC,  and 
Figures  20  through  23  give  the  new  results  for  FREQ.  The 
figures  show  that  double  exponential  smoothing  with  an  a  of 
0.80  produced  the  lowest  MAPE  for  both  SC  and  FREQ.  The  new 
low  MAPE  for  SC  of  3.17  is  lower  than  the  previous  MAPE  of 
6,92,  as  should  be  expected  after  eliminating  the  large 
errors  caused  by  periods  6  and  7.  The  new  low  MAPE  of  .95 
for  FREQ  is  not  quite  as  low  as  the  previous  of  .92. 
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Single  Exponential  Smoothing  Errors  for 
System  Constant  (Periods  1-21) 


Double  Exponential  Smoothing  Errors  for 
System  Constant  (Periods  1-21) 


Single  Moving  Average  Errors  for 
System  Constant  (Periods  1-21) 
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Double  Moving  Average  Hrrors  for 
System  Constant  (Periods  1-21) 
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Double  Exponential  Smoothing  Errors 
for  Frequency  (Periods  1-21) 


Figure  12 
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Figure  16 

Single  Exponential  Smoothing  Error 
for  System  Constant  (Periods  8-21) 
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Figure  17 

Double  Exponential  Smoothing  Errors 
for  System  Constant  (Periods  8-21) 


Figure  18 
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Figure  19 

Double  Moving  Average  Errors  for 
System  Constant  (Periods  8-21) 
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Figure  21 

Double  Exponential  Smoothing  Errors 
for  Frequency  (Periods  8-21) 


However,  there  was  very  little  confidence  in  the  previous 
results  not  only  due  to  the  two  distorted  periods,  but  also 
because  the  MAPE  was  based  on  a  single  forecast  in  period  21. 
There  was  only  one  period  forecast,  because  the  double  moving 
average  method  averaging  over  ten  periods  requires  20  periods 
of  data.  Therefore,  double  exponential  smoothing  was 
selected  as  the  best  of  the  methods  tested  to  forecast  both 
SC  and  FREQ.  The  forecasts  resulting  from  this  method,  with 
an  a  of  0.80,  are  plotted  against  the  actual  values  of  the 
indices  in  Figures  24  and  25.  These  show  that  the  forecasts 
closely  follow  the  actual  values. 

Model  Predictions 

Using  double  exponential  smoothing  with  an  a  of  0.80 
to  forecast  SC  results  in  forecasts  of  123.945  and  118.841 
for  periods  22  (March  1981)  and  23  (June  1981)  .  The  absol¬ 
ute  percent  errors  of  these  forecasts  are  10.04  for  period 
22  and  3.23  for  period  23.  The  large  error  for  period  22 
was  caused  by  the  sharp  increase  in  period  21.  The  relatively 
small  error  in  period  23  shows  that  the  forecasts  adjust 
quickly  to  large  changes  due  to  the  high  a.  The  FREQ  fore¬ 
casts  for  periods  22  and  23  are  3.776  and  3.752  with  absol¬ 
ute  percent  errors  of  0.80  and  0.79  respectively. 

These  forecasts  were  substituted  into  Models  1,  3, 
and  6,  and  VSL$  was  calculated  for  each  model  based  on  the 
8  used  in  both  periods  22  and  23.  The  model  predictions  for 
period  22  are  shown  in  Table  18  with  the  error  for  each  model. 
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Figure  24 

Actual  System  Constant  Versus  Best 
Forecast  Method  (Periods  8-21) 


TABLE  18 


Model  Predictions  for  Period  22 
(8  =  43,000) 


Model 

Actual 

VSL$ 

Predicted 

VSL$ 

Absolute 

Error 

1  Error 

1 

82.16 

88.41 

6.25 

7 . 61 

3 

82.16 

91.75 

9 . 59 

11.67 

6 

82.16 

85.44 

1.28 

1 .  56 

TABLE  19 

Mod  e  1 

Predictions  for 
(8  ■  35,000) 

Period  23 

Model 

Actual 

VSL$ 

Predicted 

VSL$ 

Absolute 

Error 

%  Error 

1 

92.59 

88.81 

3.78 

4.08 

3 

92.59 

89.79 

2.80 

3.02 

6 

92.59 

88.45 

4.14 

4.17 

The  results  for  period  23  are  given  in  Table  19.  Model  6 
has  the  lowest  error  in  period  22  and  the  highest  error  in 
period  23.  Model  3  produced  the  opposite  results,  with  the 
highest  error  in  period  22  and  the  lowest  in  period  23. 
However,  the  error  for  Model  6  in  period  23  is  still  rela¬ 
tively  low  and  close  to  the  errors  for  the  other  two  models, 
while  the  Model  3  error  in  period  22  is  very  large.  Model  1 
did  not  produce  the  best  or  the  worst  predictions  for  either 
period  and  appears  to  react  slowly  to  changes. 

The  slow  reaction  to  changes  by  Model  1  could  be  due 

93 


hi  i 


to  the  OLD  SC  term  in  the  model,  which  would  cause  the  pre¬ 
dictions  to  lag  behind  the  actual  values.  Model  3,  which 
is  dominated  by  the  SC,  appears  to  react  very  quickly  to 
changes  in  SC  and  perhaps  overreact  when  the  changes  are 
large.  Model  6  includes  the  FREQ  index,  which  seems  rela¬ 
tively  more  stable  than  the  SC.  This  appears  to  somewhat 
dampen  the  effect  of  large  changes  in  SC,  although  the  model 
still  seems  to  react  quickly  to  SC  changes.  These  results 
show  that  none  of  the  models  is  clearly  superior  in  all  cir¬ 
cumstances,  and  that  their  relative  performance  is  affected 
by  the  magnitude  of  the  index  forecast  errors. 

Summary 

Of  the  six  regression  models  tested,  three  were  re¬ 
jected  because  they  fit  the  data  poorly.  The  three  remaining 

2 

models  had  high  r  values,  and  each  was  used  to  predict  VSL$ 
for  two  periods  using  forecast  values  for  the  required 
indices.  The  results  showed  that  the  two  indices  required 
could  be  forecast  fairly  accurately,  except  when  there  is  a 
significant  change  in  an  index  trend.  The  results  of  the 
model  predictions  were  mixed  with  the  best  prediction  appar¬ 
ently  being  determined  by  the  nature  of  the  changes  in  the 
indices.  However,  with  one  exception,  the  cumulative  errors 
from  the  forecasts  and  the  regression  models  were  quite 
small . 
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CHAPTER  5 


CONCLUSIONS  AND  RECOMMENDATIONS 

Based  upon  the  results  presented  in  Chapter  4,  it 
appears  that  when  the  independent  variables  are  known,  a 
regression  model  can  define  the  VSL$  -  g  relationship  fairly 
accurately.  With  the  three  best  models,  96  to  97  percent  of 
the  variation  in  VSL$ ,  during  the  21  periods  in  the  data 
base,  could  be  explained  by  g  and  the  other  independent  vari¬ 
ables  . 

Using  one  of  these  models  will  eliminate  the  require¬ 
ment  to  calculate  several  points  on  the  VSL$  -  6  curve  directly 
from  the  current  data.  Instead,  the  points  can  be  calculated 
from  the  regression  equation. 

The  choice  of  which  regression  equation  to  use  is 
not  completely  clear  based  solely  on  the  results  presented  in 
Chapter  4,  because  these  did  not  show  one  model  to  be  sub¬ 
stantially  more  accurate  than  the  others  in  all  instances. 

Model  1,  however,  has  more  weaknesses  than  either  Models  3 
or  6.  First,  there  is  some  question  about  whether  Model  1 
actually  describes  causal  relationships  between  the  variables 
or  is  simply  the  linear  equation  which  fits  the  historical 
data  the  closest.  In  developing  this  model,  there  was  no 
underlying  hypothesis  about  the  relationships  or  even  the 
variables  to  be  included  in  the  model.  In  addition,  the  r2 
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for  Model  1  was  not  as  high  as  the  r  for  Model  6,  and  one  of 
the  coefficients  had  a  relatively  low  significance.  Finally, 
Model  1  did  not  produce  the  best  prediction  for  either  of  the 
two  test  periods.  Based  on  these  considerations,  Model  1 
appears  to  be  inferior  to  both  Models  3  and  6. 

Both  Models  3  and  6  were  derived  from  the  theoretical 
relationships  between  the  variables.  This  provides  high  con¬ 
fidence  that  the  models  define  a  causal  relationship  between 
VSL$ ,  0,  and  the  other  variables,  which  will  hold  true  for 
future  periods.  Although  either  model  might  serve  to  define 
the  VSL$  -  0  relationship.  Model  6  has  some  advantages  over 
Model  3.  One  important  difference  between  the  two  is  that 
Model  6  had  the  highest  r  .  Model  6  also  had  a  lower  signifi¬ 
cance  for  its  coefficients  than  Model  3;  however,  the  signifi¬ 
cance  levels  were  still  high  enough  to  have  a  very  high  confi¬ 
dence  that  the  variables  were  responsible  for  determining  VSL$ . 
Perhaps  the  most  important  difference  between  these  two  models 
was  in  the  prediction  errors  for  the  two  test  periods.  While 
both  models  had  the  highest  error  in  one  period  and  the  low¬ 
est  in  the  other  period,  Model  6  had  a  much  lower  average 
error  because  Model  3  greatly  overreacted  in  period  22.  This 
caused  Model  3  to  have  an  average  error  for  the  two  periods 
of  7.345  percent,  while  Model  6  had  an  average  error  of  3.015 
percent.  Both  the  higher  percent  error  of  Model  3  and  its 
apparent  tendency  to  overreact  to  changes  in  SC  were  considered 
important  drawbacks.  Therefore,  Model  6  is  recommended  as  the 
model  to  use  to  define  the  VSL$  -  0  relationship. 


The  values  of  SC  and  FREQ  which  are  needed  in  Model 
6  can  either  be  forecasted  or  obtained  directly  from  the  cur¬ 
rent  data.  The  results  showed  that  FREQ  can  be  forecasted 
with  an  average  error  of  less  than  one  percent.  The  forecasts 
for  SC  were  also  very  accurate  except  in  period  21,  when  there 
was  a  sharp  increase  in  the  actual  value.  The  error  was  also 
large  in  period  22,  while  the  forecast  model  was  adjusting  to 
the  new  SC  level.  It  is  the  opinion  of  this  writer  that  such 
large  changes  in  the  SC  level  could  be  anticipated  by  DESC  per¬ 
sonnel.  In  these  periods  when  the  SC  forecasts  will  produce 
unacceptable  errors,  the  value  of  SC  should  be  determined  dir¬ 
ectly  from  the  current  data  either  by  calculating  the  value  for 
all  items  or  using  a  representative  sample.  During  all  other 
periods,  the  double  exponential  smoothing  forecasts  should  be 
accurate  enough  to  use  in  Model  6  to  select  the  desired  level 
for  VSL$  and  8. 

The  results  of  this  research  indicate  that  the  cost  of 
various  levels  of  inventory  performance  can  be  found  through 
aggregate  measures  of  the  inventory  characteristics  using  a 
causal  regression  model  and  well-known  forecasting  techniques. 
While  the  research  concentrated  on  the  DESC  inventory  using 
the  DLA  model  to  set  safety  stock  levels,  the  results  may  be  of 
interest  to  other  DoD  components  which  also  use  the  general  DoD 
safety  level  model.  Other  organizations  responsible  for  large 
inventories  may  also  find  some  value  in  the  general  technique 
used  here  and  the  indication  from  these  results  that  the  rela¬ 
tionship  between  cost  and  performance  can  be  predicted  with  a 
fair  degree  of  accuracy. 
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APPENDIX  A 
HISTORICAL  DATA 
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Period 

Quarter 

Beginning 

0 

SC 

(mill ions) 

OLD  SC 
(millions) 

VSL$ 

(millions) 

1 

Jan  76 

50,000 

$69.60 

$76.41 

$48.86 

2 

Apr  76 

40,000 

66.74 

69.60 

46.68 

3 

Jul  76 

40,000 

65.63 

66.74 

44.62 

4 

Oct  76 

40,000 

62.92 

65.63 

43.27 

5 

Jan  7  7 

39,000 

58.03 

62.92 

40.44 

6 

Apr  7  7 

35,500 

44.26 

58.03 

33.78 

7 

Jul  77 

23,410 

39.70 

44.26 

34.30 

8 

Oct  77 

20,000 

64.90 

39.70 

47.01 

9 

Jan  78 

25,500 

66.78 

64.90 

49.88 

10 

Apr  78 

25,200 

68.82 

66. ’8 

51.30 

11 

Jul  78 

26,100 

72.57 

68.82 

55.06 

12 

Oct  78 

26,100 

78.58 

72.57 

61.57 

13 

•Ian  79 

35,000 

79.95 

78.58 

58.07 

14 

Apr  7  9 

38,000 

82.90 

79.95 

59.99 

■ 


m— mmiMw 


an 


Period 

DEMANDS 

(millions) 

ITEMS 

0- ITEMS 

FREQ 

(millions) 

1 

$256.54 

157,937 

18,733 

3.834 

2 

253.46 

156,161 

15,062 

3.826 

3 

245.78 

160,742 

17,286 

3.853 

4 

239.00 

156,596 

16,201 

3.825 

5 

214.48 

154,312 

16,710 

3.785 

6 

174.67 

105,345 

15,922 

3.147 

7 

143.51 

105,731 

14,258 

2.918 

8 

260,23 

155,153 

16,620 

3.849 

9 

270.40 

159  ,195 

15,664 

3.915 

10 

275.75 

157,562 

15,296 

3.886 

11 

288.36 

165,032 

18,296 

3.943' 

12 

309.71 

153,853 

14,714 

3.851 

13 

310.31 

152,686 

13,661 

3.803 

14 

315.20 

149,992 

17,254 

3.769 

100 


Period 

Quarter 

Beginning 

8 

SC  OLD  SC 

(millions) (millions) 

VSL$ 

(millions) 

15 

Jul 

79 

39,400 

584.82 

$82.90 

$59.07 

16 

Oct 

79 

35,000 

85.32 

84.82 

59.38 

17 

Jan 

80 

39,000 

85.00 

85.32 

65.22 

18 

Apr 

80 

46,000 

88.28 

85.00 

59.60 

19 

Jul 

80 

43,000 

90.25 

88.28 

65.00 

20 

Oct 

80 

48,200 

92.64 

90.25 

57.68 

21 

Jan 

81 

32,800 

111.63 

92.64 

84.60 

22 

Apr 

81 

43,000 

112.64 

111.63 

82.16 

23 

Jul 

81 

35,000 

115.12 

112.64 

92.59 

Period  ^ 

DEMAND$ 

millions) 

ITEMS 

0  -  ITEMS 

FREQ 

(millions) 

15 

$311.68 

150,694 

18,170 

3.705 

16 

314.14 

148,230 

15,795 

3.686 

17 

325.21 

147,695 

16,630 

3.662 

18 

305.91 

148,232 

19,376 

3.668 

19 

314.89 

151,431 

18,936 

3.687 

20 

327.50 

151,541 

20,468 

3.735 

21 

391.87 

151,720 

16,008 

3.750 

22 

390.62 

151.613 

17,156 

3.746 

23 

396.98 

150,627 

13,999 

3.724 

APPENDIX  B 

RESULTS  FROM  REGRESSING  VSL$ 
AGAINST  ln(B) 
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Period 


1 

165.00 

-29.70 

.9994 

2 

154.60 

-29.20 

.9979 

3 

152.04 

-29.08 

.9969 

4 

148.28 

-28.44 

.9975 

5 

139.60 

-27.00 

.9967 

6 

113.40 

-22.30 

.9966 

7 

90.98 

-18.20 

.9949 

8 

124.50 

-26.17 

.9940 

9 

143.09 

-28.75 

.9975 

10 

146.64 

-29.56 

.9975 

11 

162.86 

-33.04 

.9999 

12 

182.48 

-37.06 

.9999 

13 

187.53 

-36.41 

.9998 

14 

210.01 

-41.26 

.9999 

15 

212.14 

-41.65 

.9999 

16 

203.88 

-40.60 

.9997 

17 

224.53 

-43.47 

.9999 

18 

223.97 

-42.92 

.9999 

19 

233.90 

-44.90 

.9999 

20 

227.37 

-43.78 

.9999 

21 

26S.05 

-51.75 

.9997 
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